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Abstract 

This review is intended to give a pedagogical and unified view on the sub- 
ject of the statistics and scaling of physical quantities in disordered electron 
systems at very low temperatures. Quantum coherence at low temperatures 
and randomness of microscopic details can cause large fluctuations of physical 
quantities. In such mesoscopic systems a localization-delocalization transition 
can occur which forms a critical phenomenon. Accordingly, a one-parameter 
scaling theory was formulated stressing the role of conductance as the (one- 
parameter) scaling variable. The localized and delocalized phases are sepa- 
rated by a critical point determined by a critical value of conductance. How- 
ever, the notion of an order parameter was not fully clarified in this theory. 
The one-parameter scaling theory has been questioned once it was noticed 
that physical quantities are broadly distributed and that average values are 
not characteristic for the distributions. Based on presently available analytical 
and numerical results we focus here on the description of the total distribution 
functions and their flow with increasing system size. Still, one-parameter scal- 
ing theory does work in terms of typical values of the local density of states 
and of the conductance which serve as order parameter and scaling variable of 
the localization-delocalization transition, respectively. Below a certain length 
scale, £ c , related to the value of the typical conductance, local quantities are 
multifractally distributed. This multifractal behavior becomes universal on 
approaching the localization-delocalization transition with £ c playing the role 
of a correlation length. 
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I. INTRODUCTION 



This review is about mesoscopic conductors, statistics of related physical quantities and 
their scaling behavior under change of the size of the conductor. A few words about the 
terms used here shall tell the reader what to expect. 

Mesoscopic systems are intermediate between microscopic and macroscopic systems. The 
term goes back to van Kampen and was widely used after a work on "Physics of Mesoscopic 
Systems" by Imry |)9|. Mesoscopic conductors contain many (e.g. 10 19 ) elementary objects 
like electrons and atoms. On the other hand their properties cannot be obtained by using 
the thermodynamic limit (volume and particle number going to infinity with fixed particle 
density). The thermodynamic limit is a convenient mathematical device in solid state theory 
if one is interested in material constants. In particular, in macroscopic conductors the 
conductivity, a, is such a material constant and the conductance, G, is then given by Ohm's 
law, G = aW/L, where W is the cross section of the conductor and L its length. This 
Ohmic behavior is to be expected once the system size (denoted by one length L) is much 
larger than characteristic length scales (see Fig. 1). These are (i) the kinematic scale set by 
the de Broglie wavelength A or, for low temperatures, by the Fermi wavelength Xp, (ii) the 
elastic mean free path l e , which is the average distance an electron travels before its initial 
momentum relaxes and (iii) the phase coherence length L^, which is the average distance an 
individual electron travels before its initial and final phases become incoherent. Once the 
phase coherence length becomes larger than the system size the conductor can no longer 
be described by material constants. Microscopic details of the conductor will then influence 
even global physical quantities like conductance such that measurements yield finger prints 
of the mesoscopic conductor. These finger prints are due to quantum mechanical interference 
effects. 

Since mesoscopic conductors contain many electrons one has to use a quantum statistical 
description for their thermodynamic and transport quantities. However, this is not the only 
aspect of statistics in mesoscopic conductors. Apart from the electronic degrees of freedom 
the conductor contains other degrees of freedom which are responsible for resistivity to 
applied voltages. In the simplest case these are static fields caused by impurities, vacancies 
and dislocations in an otherwise ideal crystal lattice. Also static inhomogeneous electric 
fields caused by surrounding media are responsible for resistivity. One refers to these static 
fields as static disorder potentials. For different realizations of a static disorder potential the 
global conductance of an otherwise identical mesoscopic conductor will differ - sometimes 
significantly. Therefore, one considers ensembles of disorder potentials characterized by a 
few parameters, e.g. by the average mean free path on short scales. The task of the theory 
is then to determine the probability distribution of physical quantities like conductance for 
such ensembles. Thus, on top of the quantum statistical treatment for an individual disorder 
potential the statistical properties of a whole ensemble of disorder potentials is addressed in 
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mesoscopic physics. 

Typically, the phase coherence length is of the order of a few microns for metals and semi- 
conductors below liquid-helium temperatures. This is the reason that technological advances, 
starting in the 80's, were needed in order to study mesoscopic conductors systematically. 
We will discuss some of the fundamental experiments in section [TI[ 

Although it is true that technological advances were needed to bring mesoscopic physics 
to the attention of a broader community in solid state physics, a great deal of the physics had 
been discussed much earlier. For example, Anderson had introduced the idea of localization 
already in 1958 H. By localization it is meant that strong disorder can trap electrons by 
quantum interference to a finite region such that the conductor actually behaves as an insu- 
lator. Also, Landauer had presented a formula in 1957 ^54\ that describes the conductance of 
a phase coherent conductor in terms of pure quantum mechanical transmission amplitudes, 
and a similar important approach to the conductance of mesoscopic systems goes back to 
works in 1972 by Edwards and Thouless H3 . 

The works by Edwards and Thouless [|5[] and by Wegner [|131|| form the starting point 
of the scaling theory for mesoscopic conductors. Later in 1979 it was formulated in a simple 
and predictive way by Abrahams, Anderson, Licciardello and Ramakrishnan Jl[ . The scaling 
theory for disordered mesoscopic conductors is based on the assumption that the transport 
properties of systems probed at scales much larger than the mean free path should become 
insensitive to the microscopic origin of the disorder. The scaling theory states that once we 
know the conductance for a given system size we can obtain the conductance for an even 
larger system from a universal flow equation. Now, since we have already pointed out that 
the conductance is a random variable which depends in a complicated manner on the par- 
ticular disorder realization this assumption needs further explanation. As has been stressed 
by Anderson et al. U and later by Shapiro ||115yil6|| it can be correct only in a proba- 
bilistic sense. The whole probability distribution of the conductance approaches a universal 
distribution function which depends only on very few characteristic system parameters. 

The scaling theory leads to a classification of mesoscopic electron systems into three 
classes: depending on the initial values of the characteristic system parameters (called scal- 
ing variables) for a system size L the system will flow under increase of system size to an 
insulating state (localized phase), or to an conducting state (delocalized phase), or it stays 
scale independent at a critical state . Thus a disorder induced transition from a localized 
to a delocalized phase can occur which resembles a critical phenomenon. This transition is 
frequently denoted as Anderson transition or metal-insulator transition (MIT). To distin- 
guish this disorder induced MIT from other MITs occurring in solid state physics we adopt 
the notion of localization- delocalization (LD) transition. Accordingly, the critical point is 
referred to as the LD transition point. 

The central aim of this article is to provide a unique picture about the LD transition 
in mesoscopic disordered electron systems. To achieve this aim we will introduce the basic 
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physical quantities, scales and theoretical models and rely on simple calculations avoiding 
involved technical details. We stress general concepts and methods such as quantum com- 
position rules for conductors, Fokker Planck equations, and scaling laws for distribution 
functions. Important results that are based on technically more elaborate methods will only 
be quoted without detailed description. Nevertheless, we hope that the present article will 
help to find a way through the original literature. 

A number of reviews and books exist on related subjects discussed in this article. The 
following short list shall help for further reading. The review by Y. Imry |69| is an introduc- 
tion to mesoscopic physics. A review on the subject of mesoscopic transport that contains 
also information about the technological aspect of devices was written by Beenakker and 
van Houten ||15|| . The book by Datta |38] on "Electronic Transport in Mesoscopic Physics" 
covers many subjects that have proved to be relevant in the field. As review on localization 
and scaling theory for mesoscopic systems we recommend the ones by Lee and Ramakr- 
ishnan [58] and by Kramer and MacKinnon |]79| . The book by Al'tshuler and Lee J7[ on 
"Mesoscopic Phenomena in Solids" contains review articles by inventors in the field. 

The present article is organized as follows. In Sect. [TJ fundamental experiments will 
be discussed. The electronic Aharonov-Bohm-Effect makes the definition of mesoscopic 
conductors more precise. Universal conductance fluctuations demonstrate the finger-print 
character of mesoscopic conductance measurements and experiments on weak and strong 
localization point to the necessity of a critical phenomenon description of the LD transitions. 
Basic physical quantities are the linear response quantities such as conductance. These 
and the corresponding relevant energy and length scales are the subject of Sect. |T|. To 
discuss the statistical problem inherent in the problem of mesoscopic conductors we introduce 
Hamiltonian models in Sect. |V]. There we will see that the Green's function contains the 
essential information. Only the asymptotic behavior of the Green's function is needed for 
describing global conductances resulting in a scattering theoretical formulation of mesoscopic 
conductors (Sect. [V]). This point of view has been particular useful in describing strong 
localization in so-called quasi-one-dimensional systems. However, the scattering theoretical 
modeling is not restricted to this case but can be extended to more general geometries by 
considering networks of scattering matrices. After the basic quantities and tools have been 
presented we discuss the statistical properties of mesoscopic conductors in the idealized 
localized and delocalized phase, respectively (Sect. |VT|). To prepare for the statistics in 
more realistic systems we introduce the concept of the scaling theory in Sect. |V11| . The flow 
of the relevant scaling variables is controlled by so-called /5-functions from which a critical 
exponent for the LD transitions can be obtained. It describes how the correlation length 
of this critical phenomenon diverges on approaching the LD transition point. Based on a 
pedagogical model, introduced by Shapiro (cf. [|116|1 ), we emphasize that typical values of 
the conductance distribution can be good candidates for scaling variables. We pay special 
attention to the statistics of local quantities (like the local density of states) that show power 
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law scaling in non- localized phases. We come to the conclusion that, at the LD transition, 
a spectrum of power law exponents (called multifractal spectrum) is needed in order to 



describe the distribution. The more realistic models are the subject of Sect. VIII. The 



onset of broad distributions characterizing global and local quantities at large, but finite, 
values of typical conductance leaves open questions about their interpretation. To clarify 
some of them the conductance distribution in an idealized quasi-one-dimensional geometry 
is studied. The idealization is due to the assumption of absence of transversal diffusion in 
such systems. The calculations give useful insight into the physics of universal conductance 
fluctuations and lead to the interpretation that the onset of broad tails in the conductance 
distribution is intimately connected to the critical regime of the LD transition. Such regime, 
as well as long tails in the conductance distribution, is absent in the idealized quasi-one- 
dimensional conductor. Giving up the idealization, one can study the LD transition in quasi- 
one-dimensional systems by the so-called finite size scaling method. Due to the complexity 
of calculations for non-ideal quasi-one-dimensional systems the finite size scaling method is 
often restricted to numerical calculations. However, such calculations are very successful and 
provide a lot of insight into the physics as well as quantitative results for universal properties. 
The investigation of statistics and scaling of mesoscopic disordered electron systems close 
to the LD transition in Sect. IX reveals the multifractal properties of electronic states and 
opens the possibility to consider the local density of states as an appropriate order parameter 
field for the LD transition. The typical value, as given by the geometric mean, is the global 
order parameter. Similarly, we collect arguments that allow to view the typical conductance 
as the only relevant scaling variable. Our conclusions are summarized in Sect. [X] and a brief 
account of linear response theory is given in the appendix Sect. ^. 



II. MESOSCOPIC EXPERIMENTS 

In this section we will discuss some fundamental mesoscopic experiments without paying 
much attention to historical order and completeness. Rather, we wish to give motivation to 
the following sections which are based on a theoretical physicists point of view. 

The electronic Aharonov-Bohm-effect can be used for a definition of mesoscopic conduc- 
tors. It is our starting point of investigations. It points to the necessity for a quantum 
mechanical description of transport. The need to describe physical quantities in terms of 
distribution functions rather than in terms of material constants became evident after the 
discovery of universal conductance fluctuations. That one conductor can behave both, as 
a perfect conductor and as a perfect insulator, when changing system parameters is the 
insight coming from experiments on weak and strong localization. Weak localization refers 
to a situation where the system is metallic, but the conductivity is reduced as compared 
to its classical value. This behavior is interpreted as a precursor of the strong localization 
phenomenon where the system is insulating due to spatial localization of electronic states. 
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Weak localization effects have been found in various materials. The experimental verification 
of strong localization is most clear-cut in the quantum Hall effect, and first experimental 
investigations of fluctuations at the LD transition point have been performed. 



A. Electronic Aharonov-Bohm-Effect 

A ring-shaped conductor shown schematically in Fig. 2 with wire thickness of ~ 40 nm 
and a diameter of ~ 820 nm that was etched out of a high quality polycrystalline gold film 
of thickness ~ 40 nm was used to perform the first realization of the most fundamental 
experiment in mesoscopic physics, the electronic Aharonov-Bohm-Effect. At temperatures 
less than 1 K the resistance oscillates as a function of the magnetic flux enclosed by the 
ring [ 128|| . Figure 3 shows the result for the resistance oscillations found in a ring-shaped 



conductor fabricated in high mobility GaAs-Al GaAs heterostructures performed by Timp et 
al. [ |124jl . Here the diameter of the ring is 2/x m. The period in the magnetic field (~ 5 mT) 
corresponds to a periodicity in the magnetic flux enclosed by the ring of one flux quantum 
h/e where — e is the elementary charge of an electron and h is Planck's constant. The 
qualitative explanation of the oscillations relies on the concept that electrons in a magnetic 
field B = curl A, where A is the vector potential, pick up a quantum mechanical phase, 

<p=j:fA{T).dr, (1) 

c 

when they move along a classical path C that connects two points in space. Now, for an 
electron coming from the current source the path splits up into the upper and lower branch 
of the ring and the phase difference (p u — ifi between both paths is given by the total magnetic 
flux <p enclosed by the ring 

e 



h 



j> A(r)dr = ^<j> . (2) 



In terms of the the flux quantum 0o := h/e the phase difference is 27r0/0o- Taking for granted 
that the resistance is sensitive to quantum interference of electronic waves the periodicity 
of the resistance is given by 0o- The electronic Aharonov-Bohm-Effect demonstrates that 
the electrons are to be treated quantum mechanically along the ring and that their initial 
phase is not randomized before leaving the ring. This phenomenon is called phase coherence 
and the Aharonov-Bohm-Effect gives rise to an operative definition of the phase coherence 
length (see also Sect. |111B| ). On increasing the temperature or the diameter of the 



ring the phase coherence gets lost and the Aharanov-Bohm effect disappears. The system 
size at which the Aharanov-Bohm effect (gradually) disappears quantifies the temperature 
dependent phase coherence length L^. The disappearance of the Aharanov-Bohm effect 
is caused by an increasing number of inelastic scattering events which destroy the phase 
coherence. The average distance of inelastic scattering events can be identified with L^. 
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B. Universal Conductance Fluctuations 



Measurements of conductance in mesoscopic systems that have been performed since the 
80's, showed sample specific, reproducible, statistical fluctuations when system parameters 
like Fermi energy, magnetic field or disorder configurations were changed ||130|| . In Figure 4 



the results of a recent experiment are shown. Gold-nanowires of very high purity with 
cross-section 30 2 nm 2 and lengths between 400 and 1000 nm were used. The conductance 
G was measured as a function of applied magnetic field below 16 T in steps of 1 mT at 
60 niK. Although the average conductance varies from « 3000 e 2 /h (for L = 400 nm) to 
~ 1400e 2 //i (for L = 1000 nm), the variance of the conductance stayed constant and was 
~ 0.1(e 2 /h) 2 . This phenomenon is referred to as universal conductance fluctuations (UCF). 
In a mesoscopic conductor with high conductance (in atomic units e 2 /h) the conductance 
fluctuations are reproducible finger-prints of the conductor. However, the variance turns 
out to be approximately a constant 

var(G) = const. ~ 0(1) x {e 2 jhf . (3) 

Two aspects of UCF are interesting: reproducibility and universality. That the fluctuations 
are reproducible can be understood again by a semiclassical description of conducting elec- 
trons. Consider electrons moving along classical paths between elastic collisions, but having a 
quantum mechanical phase attached. Since the phase information is not lost over distances 
of the phase coherence length (mesoscopic regime) a sample specific interference pattern 
arises. The interference pattern will be altered e.g. by applying a magnetic field. This 
gives rise to reproducible conductance fluctuations. The fact that the order of magnitude 
of these fluctuations turns out to be independent of the average conductance and does only 
depend weakly on dimensionality and symmetry properties of the conductor is, however, not 
easy to understand without having a transport theory of mesoscopic conductors. A striking 
consequence of the phenomenon is that mesoscopic conductors can not be characterized by 
an average conductivity being a material constant, even if the system size L < becomes 
very large. Conductivity a is related to conductance by Ohm's law. For a cubic geometry 
the conductance behaves as G = o~L d ~ 2 . For UCF the relative fluctuations of conductance 
are SG/G oc L 2 ~ d which is in contrast to the classical behavior SG/G oc L~ d l 2 . Thus, even 
in 3D metallic systems the relative fluctuations are much stronger than 1 / \j volume. 



C. Weak And Strong Localization 



Early measurements of conductivity of amorphous Silicon |19| showed that, at low enough 
temperatures, the temperature behavior of conductivity follows Mott's T~ 1//4 -law, 



cr(T) = o- exp (-(Tq/T) 1 /^ 1 )) for d = 3 



(4) 
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over two orders of magnitude with some constants To and <To (see Fig. 5). This behavior of the 
conductivity could be understood when assuming that the transport is mediated by phonon 
assisted hopping processes between localized states |[L00|| . Localized states are spatially 
localized to some finite volume within the conductor. The finite volumes are characterized 
by a localization radius £ r . The phonon assisted hopping is possible because localized states 
can be energetically close to each other if the corresponding localization centers are separated 
by distances much larger than the localization radius. The observation of Mott's law made 
it clear that electronic states can indeed be localized by disorder which has been pointed 
out by Anderson || much earlier. The effect is referred to as strong localization since the 
localization radius is smaller than system size. 

A further experimental hint to localization came from weakly disordered metallic films at 
low temperatures. They showed a logarithmic increase of resistance when the temperature 
is decreased (see Fig. 6). Although the system behaves metallic with large conductance 
the conductance was lower than its classical value. The effect was interpreted as being 
due to quantum interference corrections to classical transport and is referred to as weak 
localization. In weak localization the localization length is larger than system size but the 
effect is interpreted as a precursor of strong localization, i.e. under further increase of system 
size (while still being mesoscopic, i.e L < L^) the system would show strong localization. 
This is, of course, in practice very hard to realize when localization lengths are larger than 
microns. 

The weak localization effect (predicted by Gor'kov et al. can be explained by the 
phenomenon of enhanced backs cattering due to quantum interference. It relies on the fact 
that the quantum coherent probability of the sum of two probability amplitudes \A + A\ 2 = 
4\A\ 2 is two times larger than the incoherent sum of the corresponding probabilities, \A\ 2 + 
\A\ 2 = 2\A\ 2 . Consider the chance for an electron traveling in a phase coherent conductor 
(see Fig. 1) to return back to some point within the conductor. In Fig. 7 a possible path is 
drawn for such a path (straight line). If the motion is time reversal symmetric also the time 
reversed path (broken line) is possible. The amplitudes of both paths add coherently in a 
semiclassical calculation of the return probability. Thus, the quantum mechanical return 
probability is twice the classical value. Without explicitly calculating the conductivity it 
is obvious from this observation that the conductivity must be reduced as compared to its 
classical value. 

The coherence is limited by the phase coherence length which varies as some inverse 
power of temperature, L^(T) oc T~ p For system sizes L larger than the phase coherence 
length L^iT) the logarithmic corrections to the classical conductivity could be interpreted 
as being due to the quantity — lnL^// e where l e is the average elastic mean free path. 

Magnetic fields destroy the time reversal symmetry of the electronic motion and thus 
the weak localization effect vanishes (see upper set of Fig. 8) A very interesting effect is due 
to spin-orbit scattering processes that are time reversal symmetric. In this case the time 
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reversed path, however, gives rise to destructive interference in the return probability. As 
a result the sign of the conductance corrections is reversed (see lower set of Fig. 8) and the 
enhancement of conductivity above its classical value is denoted as weak anti- localization 
(first predicted by Hikami et al. [ 33 1 ) . 

Weak localization experiments are an important tool e.g. to determine the phase coher- 
ence length and the spin-orbit scattering length. For more details we refer to the review 
by Bergmann |18] . 



To reach experimentally the regime of strong localization where the localization radius is 
less than the system size L < one has to use strongly disordered systems or work in energy 
regions with very low density of states. Localization-delocalization (LD) transitions have 
been observed in several experiments (for reviews see |104j]113|| ). Figure 9 shows the finding of 
an experiment by Stupp et al. ||122|| indicating that the conductivity vanishes on approaching 
a critical charge carrier concentration. A major problem in comparing theoretical models 
with the experiments comes from the fact that in the low density of states regimes Coulomb 
interaction between electrons becomes essential, but most theoretical models are based on 
an independent electron picture. 

The most clear-cut experiment for the occurrence of LD transitions is the quantum 



Hall effect discovered by von Klitzing ||76|| . The quantum Hall effect that occurs in two- 
dimensional electron gases in the presence of strong perpendicular magnetic fields is char- 
acterized by a step-function like behavior of the Hall conductivity ou as a function of the 
so called filling factor (a dimensionless quantity proportional to carrier concentration and 
inverse magnetic field) and by a vanishing dissipative conductivity a in the Hall plateau 
regimes (see Fig. 10; Note, that for non- vanishing Hall conductivity a vanishing dissipative 
conductivity causes also a vanishing longitudinal resistivity). The peaks in the dissipative 
conductivity as a function of the filling factor have a clear interpretation in terms of LD 
transitions. However, no metallic phase exists. All states, except for those at critical fillings 
where the localization radius is larger than system size, are strongly localized. The sys- 
tems showing the quantum Hall effect (quantum Hall systems) are well suited to study the 
properties of critical states. Recently, in an experiment by Cobden and Kogan it was 
possible to extract the whole conductance distribution at criticality for a truly mesoscopic 
quantum Hall system. The main finding was a conductance distribution which is indepen- 
dent of system size (within the mesoscopic regime), almost uniform between and e 2 /h. 
This means in particular that the fluctuations are of the same order as the average value. 
Thus, fluctuations in the conductance can become much stronger than one could expect 
from classical transport theory where the relative fluctuations scale like the square root of 
the inverse volume. 
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III. BASIC PHYSICAL QUANTITIES 



This section prepares for the discussion of transport theory in mesoscopic electron sys- 
tems. We introduce the basic physical quantities of disordered mesoscopic systems for which 
the phase coherence length is larger than the system size. They are introduced here on a 
phenomenological level as response quantities to applied fields. The modeling of mesoscopic 
systems and the tools to calculate these quantities is the subject of the following Sect. [IVJ 
(and the appendix The basic physical quantities are the thermodynamic density of states 
and the transport quantities resistance, conductance and diffusion. Related physical scales 
are thermodynamic and kinematic scales like level spacing, Fermi energy, Fermi wavelength, 
magnetic length and transport scales like Thouless energy, elastic mean free path and the 
localization length. 



A. Linear Response Quantities 

The thermodynamic density of states 

M - 1 (5) 

determines the change of the global particle density n with respect to a change in the 
chemical potential fi. The local thermodynamic density of states p(r; /x) is defined similarly 
when replacing n by the local particle density n(r). 

In charge transport one considers a stationary though non-equilibrium situation. Global 
transport quantities are the resistances in a multi-probe setup 

7tkl,mn = J , I = h = ~h , I&{k,l} = (6) 

where probes are characterized by electro-chemical potentials t/j. The current / is driven 
from the current source (probe k) to the current sink contact (probe I) (cf. Fig. 11). 
The resistances IZm mn are determined by the conductance coefficients G^n PI 



h = Y, G UUk-u n ). (7) 

n 

With the condition of total current conservation J2k h = the conductance coefficients Gk n 
fulfill the relations Q 



1 Equations (|6|,0) can also be considered for time dependent transport with the quantities U m , 
Ik being the Fourier components for frequency uj. In that case TZki,mn and G^m will depend on u, 
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^ Gkn — E Gkn — . (8) 

k n 

In a two-probe geometry the transport is determined by the two-probe conductance G 
(Ohm's law) 

G = I/U = R- 1 (9) 

where U — U\ — U% and G = Gi 2 and R = TZn,i2- 

A full description of the local charge transport is possible with the help of the conductivity 



tensor, a M!/ (r, r '; u), 



3 A 



J dV^^(r,r»^(r') (10) 



relating local electric field, E(r), to the local current density, j(r). Both can here be con- 
sidered as the Fourier components for frequency u of the corresponding time-dependent 
quantities (~ e~ lujt ). 

Local currents can be also generated by charge density gradients under conditions of local 
equilibrium, i.e. the chemical potential /i is a function of position. The relevant transport 
coefficient in this case is the diffusion function, D nv {r, r';u>), 

Mr) = -E/ d d r'D, u (r,r';u)^jq(r') (11) 

This name is justified since the continuity equation between current and charge density, 
q(r), div j(r) = iuq(r) tells that g(r) obeys a non-local diffusion equation 

Under the conditions of local equilibrium the charge density can be related to the gradient of 
the chemical potential via the thermodynamic charge response, n(r, r') = e 2 5n(r)/5fi(r '), 



-icoq(r) =Y.j dd ^ D ^>)^') (12) 



9 I \ I \-\ f M ITT/ l\ 9 

-q(r) — ' 1 ' ' " 1 "- 1 



dr ^-, v e)- 1 Jd d r'Tl{vy)^{v'). (13) 

Since the linear response quantities conductivity and diffusion belong to the same system 
and describe currents in response to local electric and chemical potential gradients, respec- 
tively, they can be related by identifying the currents and driving fields. This constitutes a 
generalization of the important Einstein relation between conductivity and diffusivity, 



too. The total current conservation and Eq.(|8|) will still be valid as long as electric fields which are 
generated when applying voltages to the system (including the coherent probe region) have their 
sources and sinks within the system (see e.g. pl|). 
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<V(r,r';u;) = J d d r"D^{r,r";u)U{r",r') (14) 
Equation ([K]) reduces to Ohm's law for spatially averaged currents and electric fields 



and for the global conductivity tensor, a^ u , Eq. (|13|) reduces to the global Einstein relation 
between conductivity and diffusion 

= e 2 D flu (uj)— . (16) 

In the following we will concentrate on the d.c. limit of transport (uj — > 0) and on a 
longitudinal (dissipative) component of the conductivity tensor, say o := o xx . In the case 
where off-diagonal components of conductivity vanish, e.g. in the absence of magnetic field, 
the two-probe conductance G of a system with box geometry (characterized by length L and 
cross section Lf 1 , see Fig. 12) is related to the conductivity a by 

G = aL d t ~ 1 /L. (17) 

Equation (|TT|) is also referred to as Ohm's law. Whenever off-diagonal components of con- 
ductivity are present Eq. ( |TTD refers to the dissipative conductance which can be measured 
in an appropriate four-probe setup. We stress that the global conductivity is treated here 
as a spatial average over the conductivity tensor, o~ = f d d r J d d r'a(r,r')/(LL d ~ 1 ), and may 
depend on L, L t . 



B. Relevant Scales 

The Einstein relation, Eq. ([TBI) , and Ohm's law, Eq. (|T7|) allow for a very interesting 
interpretation of the conductance measured in atomic units e 2 /h. Defining a diffusion time 
to by tp := L 2 / D and a quantum mechanically related energy scale 

E Th :=h/t D , (18) 

called Thouless energy, and the thermodynamic level spacing, A(/i), by 

(LLt 1 )^ =: (A(^))- 1 (19) 



we obtain the Thouless formula [123] for the dissipative conductance, G = (e 2 /h)g, 

g = E Th /A . (20) 
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This means that the dimensionless conductance, g, can be expressed as the ratio of two 
characteristic energy scales of the problem: a transport related energy scale and a 
thermodynamically defined energy scale A. As we will see later this expression will help a 
lot in classifying mesoscopic electron systems. 

So far we have not introduced any particular dynamical properties of the electron system. 
In a linear response treatment as outlined in the appendix Sect. [A] thermodynamic as well as 
dynamic response quantities can be calculated by appropriate correlation functions. In a first 
attempt we make a phenomenological ansatz to describe the global dissipative conductivity 
u thereby making contact with standard macroscopic transport theory. 

In the absence of magnetic fields and in the presence of static disorder one can make 
the plausible assumption that an electron starting with velocity v will decelerate owing 
to scattering by impurities. This process can be modeled in the time dependent velocity 
correlator (see Eq. ( JA4] , |A8| ) 



(v u ;v^t)) = v 2 Q e~ t/r 5 flu (21) 

where r is a phenomenological momentum relaxation time and v 2 . the static velocity auto- 
correlator for a given direction and serves also as a measure of the kinetic energy of the 
initial state. With the help of the Kubo relation, Eq. ( |A11| ), this correlator is found to be 



•* = £ (22) 

(N being the number of electrons, (3 the inverse temperature 1/ksT and m being the elec- 
trons mass). The global conductivity is proportional to the Laplace transform of the velocity 
correlator (see Eq. ( |A19p ) which yields the Drude formula for the d.c. conductivity 

e 2 nr 

&xx = , a yx = 0. (23) 

m 

For non-interacting electrons, at T = 0, Eq. (|23"D can be interpreted as 

a„(Ep) = e 2 n(E F )r(E F )/m = e 2 p(E F ) ^^ (24) 



where E F is the Fermi energy, vp = J(2/m)E F the Fermi velocity and p(E F ) the density of 
states. 

In a finite system of volume L d the density of states defines the level spacing 

A = (pL*)- 1 (25) 

which is the smallest quantum kinematic energy scale of the system. 

The Fermi velocity sets a quantum kinematic length scale, the Fermi wavelength 

X F = — (26) 
Vpm 
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The momentum relaxation time is equivalent to an elastic mean free path, l e , that describes 
a typical distance between elastic collisions 



h = V F T . (27) 

The corresponding diffusion constant follows from the Einstein relation, Eq.(|l6|). 

D= V -^. (28) 

Assuming that in a large but finite system the diffusion constant is indeed a constant the 
Thouless energy, Eq. flT8|), and the conductance, Eq. (|T9D, can be determined as 

„ hD S d ( L \ d - 2 l e 

Here Sd is a geometrical number, the surface area of the unit sphere in d dimensions, e.g. 
2n in 2D. The situation with l e ^> \ F corresponds to large g and will be denoted as weak 
scattering case. 

Within macroscopic transport theory the transport is determined by material constants, 
e.g. the diffusion constant, and the basic task is to calculate them from microscopic models. 
The situation changes when the system becomes smaller than the phase coherence length 

Consider a given conducting system with fixed geometry and at fixed temperature. Imag- 
ine that we can pierce a hole in the conductor and put a magnetic flux through it. As 
quantum mechanics tells, such flux can only influence the phase of the electrons wave func- 
tions. In case that one observes periodic oscillations in the resistances when changing the 
flux over several flux quanta, h/e, (electronic Aharonov-Bohm effect, see Sect . [I I A|) the elec- 



tronic system is mesoscopic, i.e. the electronic wave functions are phase coherent over the 
entire conducting system. As soon as the charge carrying wave packets experience scattering 
processes with dynamical degrees of freedom, e.g. phonons, that can dissipate energy from 
the electronic system the phase coherence will be destroyed and the electronic Aharonov- 
Bohm effect gets lost. It is worth noticing that elastic scattering cannot destroy the phase 
coherence. Thus, the criterion for an electron system to be mesoscopic is set by the phase 
coherence length being larger than the system size L, 

L^>L. (30) 

Here L$ can be considered as a typical distance before the charge carrying states experience 
an inelastic scattering event. It is evident that will usually decrease with increasing 
temperature. In practice one has to go below liquid Helium temperatures and use micro- 
structured electron systems to reach the mesoscopic regime. 

The main consequence of the phase coherence in mesoscopic conductors is the fact that 
quantum interference plays a dominant role for physical quantities to be observed. As 
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we have already seen in Sec. |H| ft can cause two important effects: unexpectedly large 
fluctuations and localization. 

The first effect can be qualitatively described as follows. A wave packet entering the 
mesoscopic conductor will be scattered by some elastic collisions. After a while it will form 
a complicated interference pattern. If one changes an external parameter only slightly, e.g. 
magnetic field, one has to be aware of drastic effects on the measured physical quantity since 
the interference pattern changes. A similar effect can be caused by a slight modification in 
the arrangement of elastic scatterers. At least, one cannot expect that a measured physical 
quantity is close to an average value obtained from a large number of measurements under 
slightly different conditions. Instead, the quantities will show fluctuations which can be- 
come unexpectedly large. However, such fluctuations are reproducible. The magnitude of 
mesoscopic fluctuation effects will depend on system specific parameters. This will be the 
topic of Sees. |VlJjjVllJ| , p[X| . What can be learned by this reasoning is that a characterization 
of a mesoscopic sample is incomplete unless the whole distribution function of the physical 
quantity of interest is obtained. It is no longer possible to describe it by material constants 
(see Sect. |I B|) . This insight forms the motivation to ask for the statistical properties of 
mesoscopic electron systems. 

The second effect, localization, was pointed out by Anderson in a by now famous paper 
||. He demonstrated in a model calculation that the conductance of a phase coherent disor- 
dered (mesoscopic) system can decrease exponentially with increasing system size. Denoting 
L as the linear system size, the conductance can scale as 

g(L) ex exp(-2L/0 (31) 

where £ is called the localization lengttft. The origin of this phenomenon lies in the multiple 
interference of the electrons wave packet after several elastic collisions that can lead to 
destructive interference such that the wave packet is bound to a finite region, the localization 
volume, with exponentially decreasing amplitude outside of this region. For system sizes 
that are several times larger than the localization length the conductor actually behaves 
as an insulator. The effect becomes even more interesting since it is possible that one 
system can change from the insulating, localized, phase to a conducting, delocalized, phase 
when parameters are changed continuously. Such parameters can be the degree of disorder, 
electron density or external fields like pressure, electric or magnetic field (see Sect. [II Q ). 
Thus, mesoscopic systems can undergo a disorder induced metal-insulator transition which 
we called localization-delocalization transition (LD transition). 



2 The factor of 2 is due to conventions which will be discussed later. 
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IV. MODELING OF MESOSCOPIC SYSTEMS 



In this section we introduce model systems appropriate for disordered mesoscopic sys- 
tems and describe how to obtain relevant physical quantities. Since we do not want to 
describe coupling of electrons to dynamic degrees of freedom that can dissipate energy from 
the electronic system we ignore from the very beginning electron-phonon interactions and 
restrict ourself to zero temperature, unless otherwise stated. In addition, we will treat elec- 
trons mainly as independent additive systems characterized by effective mass, only. This 
is a drastic restriction and it is in general not justified to incorporate electron-electron in- 
teractions into such effective one-particle models. However, with one-particle models one 
can already gain some insight into the mechanism of LD transitions and fluctuation effects. 
The interplay of disorder and electron-electron interaction forms a challenging research field 
which has still to be developed. Furthermore, internal degrees of freedom of the electrons 
will be ignored to a large extent, leading to models of non-interacting spin-less fermions in 
the presence of static disorder potentials. 



Independent electrons of effective mass m moving in a random one-particle potential 
V(r) are described by the Hamiltonian 



where the kinetic momentum p + eA may include magnetic fields B = curl A. Randomness 
of the disorder potential has to be specified by a distribution function. It is often assumed 
that the essential parameters of the distribution are contained in the first and second cu- 
mulant, i.e. a Gaussian distribution, fixed by the mean value (V(r)) and a finite range 
correlation function (V(r)V(r')), will be general enough to cover the essential physics. It is 
reasonable to assume that neither specific points, nor specific directions are preferred by the 
whole ensemble of potentials. Then one can characterize the ensemble by two parameters, 
a potential correlation strength, Vq, and a potential correlation length, ly, 



In case that ly can be considered to be smaller than other relevant microscopic scales, it 
can be set to zero, leading to a Gaussian white-noise potential 



A. Hamiltonian Modeling 




(32) 



(V{r)) = 
(V(r)V(r')) =\/ 2 exp(-|r-r'| 2 /^) 



(33) 
(34) 



(V(r)V(r')) = 7 2 5(r-r') 



(35) 



18 



where 7 2 corresponds to {^/tiIvYVq . Within the Gaussian white-noise potential model the 
disorder is characterized by one parameter, 7 2 , of dimension [Energy 2 Volume] . We will see 
later, Eq. (|59f) , how this parameter corresponds to an elastic mean free path, l e . 

A homogeneous magnetic field pointing in a fixed direction, say B = (0, 0, B), introduces 
a quantum mechanical length scale, the magnetic length l B defined by the size of an area 
which is penetrated by one flux quantum h/e 



l B = ^h/{eB) . (36) 

The magnetic length comes along with chirality due to the axial vector character of magnetic 
field. As long as the magnetic length is larger than the actual system size, the effect of the 
magnetic field on the classical paths of electrons can be neglected and it will only influence 
the phases of wave functions. As soon as the magnetic length is of the order of the mean 
free path l B ~ l e the dynamics is, to a large extent, determined by the magnetic field and 
we refer to this situation as strong magnetic field case. 

A simple modeling to account for possible spin-orbit scattering events that can flip the 
electron's spin is given by the addition of 

S-(gradf/ so xp) (37) 

to the Hamiltonian in Eq. K5,|52|. Here S is the spin operator and Uso denotes an 
average spin-flip potential the gradient of which is taken to be a constant into all directions 

grad£/so = {au,vu, °u) ■ (38) 

Depending on the kinematic properties of the electron's state the inverse of crjj gives rise 
to a distance, the spin-orbit scattering length Igo, after which the electrons spin relaxes on 
average. 

So far, we have considered a Hamiltonian modeling in a continuous real space represen- 
tation. The randomness is characterized completely by a diagonal potential, Eqs. Q33| , |53] , |55D . 



It is instructive to consider the same model in a discrete matrix representation and to ask 
for the resulting statistical properties. In a finite cubic system (with volume L d ) with peri- 
odic boundary conditions we can use the momentum representation, p | k) = hk | k), which 
diagonalizes the kinetic energy in Eq. (|3"2|)(B = 0) 



tf kk , = ^<W+V(k-k') (39) 

where V(k) is the Fourier transform of V(r). The Gaussian white noise results in the 
correlation 

V-(k! - k 2 )V(k 3 - k 4 )\ = L-V^-ka.ka-k. • ( 4 0) 
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Thus, the random matrix has a diagonal deterministic part and a random part the 
average value of which vanishes, but with strong statistical correlations among those elements 
which fulfill momentum conservation. What can be seen here is that the statistical properties 
of a representing matrix for the Hamiltonian do, in general, depend on the chosen basis in 
Hilbert space. 

The matrix can be taken to be finite dimensional if one restricts the wave lengths 
to be larger than a convenient microscopic scale, a. Then the matrix dimension is ~ (L/a) d 
and the system can be studied on a computer. 

A direct modeling of disordered mesoscopic systems by a finite dimensional Hamiltonian 
matrix is given by a tight-binding version of the Gaussian white-noise potential and was 
introduced in the work of Anderson ||. This Anderson model is defined on a cubic lattice 
with lattice constant a 

# = ^£ m |m)(m|+ ^m,m'| m )( m 'h ( 41 ) 

m <m,m') 

where |m) denote tight-binding states situated at lattice point m and (m, m') means nearest 
neighbors only. The site energies e m correspond to the (random) potential energy and the 
(deterministic) hopping matrix elements £ mim ' to the kinetic energy. Again, the matrix 
dimension is ~ (L/a) d for a cubic system of linear size L. However, the matrix is sparse 
since the nearest neighbor condition leaves most elements vanishing. 

In the absence of magnetic field the hopping strength is taken to be equal for all nearest 
neighbors (with coordination number Z) and set to 1. This defines the kinetic energy band 
to be of width 2Z. 

A homogeneous magnetic field B = (0, 0, B) can be included in the kinetic energy by the 
Peierls substitution | |106| |. In the Landau gauge, A = (0, — Bx, 0), it reads in 2D 



tm,m'(B) — t 

m,m' 

(0) {e ±2 — 

Om' x ,Tn x Vm' y ,m y ±l } • (42) 

Here < a < 1 is the number of flux quanta h/e per unit cell. In contrast to continuum 
models a feature of commensurability appears. Requiring commensurability between the 
lattice constant and the magnetic length restricts a to take rational values only 

a = p/q . (43) 

Furthermore, in this model the magnetic length is restricted to be larger than the lattice 
constant since the effect of magnetic field is assumed to be weak on the size of a unit lattice 
cell, i.e. only its influence on phases is taken into account by the Peierls substitution. 

The statistical properties of the model are fixed once a distribution of the independent 
diagonal energies is chosen. A convenient choice is a box distribution on the interval [— Vo, Vo] 
with correlation 

V 2 

{e m e m >) = y^mm' (44) 
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The parameter V$ /3 can be identified with 7 2 / (^/waY (7 taken from the white-noise model 
of Eq. (^)) and determines the mean free path l e (see Eq. (^)). We mention that the 
inclusion of spin-orbit scattering into the Anderson model is also possible |48| , p!o|j5T[| and 
gives rise to a spin-orbit scattering length l$o- 



Equations (|3~2|j41"D correspond to the most prominent models of disordered mesoscopic 
electron systems. In the following we want to discuss, how one can use these models to 
calculate some of the basic physical quantities. 



B. Physical Quantities by Green's Function 

As outlined in the appendix, for non-interacting Fermion models we can concentrate on 
zero temperature (see Eq. ( |A27| )) and only need the resolvent of the Hamiltonian 

G ± (E) = (E ± - H)-\ E ± = E±ie, e -f +0 (45) 



to calculate the basic physical quantities (Eqs. (|A24| - |A32|) ). Choosing a convenient basis 



one has to perform traces involving G ± {E). The knowledge of G ± {E) is equivalent to finding 
the eigenvalues e a and eigenvectors | ip a ) of H since 

The dynamic content of the resolvent is expressed by its relation to the time evolution 
operator U(t) = e ~' lHt l h 

±00 

ihG ± {E)= J dte {iETe)t/h U{t) . (47) 


With this it can be concluded that the long time average of the time dependent probability, 
P(r, r '; t), to find an electron at point r ' that was created at point r is related to the Green's 
function, G ± (r,r';E) = (r' | G ± (E) | r), by 

P(r,r';t)= hm^- J dE \ G + (r,r';£) | 2 . (48) 

That leads to the important interpretation that (e/27r) | G + (r, r E) | 2 defines the long time 
transition probability from r to r ' for a given energy E. 

The local density of states (LDOS) at T = 0, p(r; E F ) = (r | 6(E F — H) \ r) can be 
calculated by (see Eq. ( |A22j )) 

p(r; Ep) = -vr^Im G + (r, r; E F ) = - f d d r' \ G+(r, r E F ) \ 2 = ^ | V«(r) | 2 S(E F - e a ) 

(49) 



21 



and the global (averaged over volume L d ) density of states (DOS) by 

p(E F ) = L- d Y.S(E F -e a ). (50) 

a 

For a finite closed system the DOS consists of isolated 5 peaks which have to be smeared 
out by attaching a finite width of the scale of the level spacing to the peaks or by averaging 
over the ensemble of disorder realizations. The average level spacing A at E-p can thus be 
obtained by 

A = -7r(lm TrG+V 1 (51) 



Furthermore, it is of conceptual importance to mention that the Green's function 
G(r, r '; E) can be represented by a path integral, since 

(r ' | U(t) | r) = J d[r(t)} exp(i/h)S[r(t)] , (52) 

where the integration runs over all paths r(t) connecting initial (r) and final (r') points 
and 5*[r(t)] denotes the classical action associated to the Hamiltonian of the electron. The 
stationary path of 5'[r(t)] yields the solution of the classical equation of motion. 

As a first instructive example we calculate the ensemble average of the Green's func- 
tion in the continuous model without magnetic field, using the momentum representation, 
Eqs. (^,^). Since we cannot solve the Schrodinger equation for arbitrary disorder potential 
we rely on a perturbative analysis, 



oo 

G+ = G +°J2(VG +0 ) n , (53) 



n=0 



where (k'|G ±0 (£) |k) = (E± - {Tikf / 2m)~ l 5^, is the Green's function of the kinetic energy 
term. Due to the property Eq. ([40]) the ensemble averaged Green's function is diagonal in 
momentum representation, too. It can be written in the form 



(k| G ± (E) |k) = (E ± - (hkY/2m - S ± (k)J (54) 

where S ± (k) is called the self-energy the imaginary part of which defines an inverse time 
scale which can be interpreted as an elastic mean free time r 

|- = ImS + (k). (55) 

We mention that, in general, this r must not be identical to the average momentum relax- 
ation time due to random collisions. 

The self-energy can now be calculated order by order in perturbation theory (for details 
cf. [0]). For weak scattering, i.e. for l e ^> X F , it turns out to be sufficient to take the first 
non-vanishing contribution which yields 
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E + (k) = rVE(k|G + °|k) (56) 

k 

and improve it by replacing G +0 with the full resolvent G + in Eq. ([56]) making the approx- 
imation for the self energy self-consistent. It leads to 

R (57) 



2-K r y 2 v 



where v is the average DOS which, for Fermi energies far off the kinetic energy band edges, 
is not changed as compared to the DOS without disorder, i.e. 

mS d (2n\ d - 2 
U{Ef) = (2^F (58) 

Consequently, the elastic mean free path is related to 7 and by 

l e = v F r = — - oc — - — Xf ■ (59) 
m 7 iJd 7 

Calculating the disorder average of the conductivity requires already the average over 
a product of two Green's functions (called two-particle Green's function) which is already 
a non-trivial task. Starting from the Kubo formula, Eq. ( |A28|) , for the velocity correlator, 
using the momentum representation and replacing the matrix elements of the resolvent 
by their average values leads back to the Drude formula with the momentum relaxation 
time given by Eq. Q5"7|). However, this procedure ignores completely that the average of 
two Green's functions does not decouple into a product of averages. In the context of the 
path integral representation for the Green's functions, Eq. fl52"D, it becomes obvious that 
this method cannot cope with quantum interference effects. Thus, the mean free time as 
given by Eq. Q57| ) corresponds to a Drude like conductivity when quantum interference is 
suppressed. 

The symmetry properties of the disorder ensemble yield a translational and rotational 
invariant two-particle Green's function 



2e 2 e 2 

a - 



|G+(r,r',£)| 2 ) =: K(\r — r'\,E) . (60) 
With this function the longitudinal conductivity, as given by Eq. ( |A32j ), is 

J d d rx 2 K(\r\,E F ) . (61) 

The diffusion constant D(E) is defined by the long time limit of the square displacement 
with respect to the probability distribution eK(r, E) as 

p(E)D(E)t(e) := e J d d r x 2 K(\r\, E) (62) 
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where t(e) = h/2e is a growing time scale as e is sent to +0. Thus, Eq. (|6lD demonstrates 
the validity of the Einstein relation. 

Furthermore, it allows for a first definition of a localization length by means of the 
Green's function, i.e. if eK(\r\,E) stays finite for arbitrary finite r, but falls off as 

etf(|r|,£)ocexp(-2|r|/f ) (63) 

the average conductance of a large volume L d will show the localization phenomenon intro- 
duced in Eq. (|31|) with £ defining a localization length. Note, that this length is defined 
with respect to the average of \G + \ 2 and must not be identical to just 1/q of the exponential 
fall off of averages of powers |G + | 2|J . Nevertheless, the factor of 2 in Eq. (^) follows the 
convention to associate the localization length with the modulus of the Green's function. 



C. The Statistical Problem 

We are now in the position to formulate the statistical problem of calculating physical 
quantities in mesoscopic electron systems. The system is described by an ensemble of Hamil- 
tonian matrices in a certain matrix representation. The relevant physical quantities can be 
obtained from the Green's function G + (r,r';E). The average value will only determine 
the average DOS and the short distance mean free path, however not the global transport 
behavior. To account for the latter, the averaged two-particle Green's function is needed in 
order to see localization effects. Fluctuation effects in global transport cannot be obtained 
on the basis of averaged two-particle Green's functions; at least four-particle Green's func- 
tions have to be considered. We see that this situation calls for a general approach to the 
distribution functions of G + (r, r '; E). The field theoretic approaches to mesoscopic systems 



commonly known as non-linear sigma-models are constructed to fulfill this purpose ||132| , f46 
We will not enter this subject because of its technical complexity, but will here focus on 
some general considerations. 



Owing to Eq. ([46]) the statistics of G + (r,r';E) is contained in the joint probability 



distribution of eigenvalues, e a , and eigenvectors, ip a , 

i Ipai ? £«2 ? ? ^-az i ^03 ) • • •) ■ (64) 

Of course, it seems hopeless to find a general method to determine this function for any 
given ensemble of Hamiltonian matrices. 

To get some feeling about the nature of the statistical problem we consider the Hamil- 
tonian in a finite basis {|«)} i=1 N as hermitean N x N matrix = (i | H | k). The 
diagonalizing unitary matrix U G U(N) with matrix elements U^a fulfills 

H UjuHikUka = e a S a(3 . (65) 
ik 
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It is related to the amplitude of an eigenstate ipa in the {\k)} basis by 

Mk) ■= (A# Q > = U ka (66) 

We can think of the randomness of H as being controlled by a large number of independent 
parameters, e.g. strength and position of point scatterers in space, or the on-site energies in 
the Anderson model. Each realization of the Hamiltonian represents one point in this high 
dimensional parameter space. We can ask what happens for a certain stochastic process 
H(s) in this parameter space where s denotes a fictitious time and each element H(s) is 
weighted by the corresponding probability. The original H corresponds to some arbitrary 
point, say H(sq). We can increase s in small steps 5s which leads to 

H (s + 5s) = H(s) + 6H(s, 5s) (67) 

The corresponding unitary matrix U (s) behaves then as 

U(s + 5s) = U(s)U(s,5s) (68) 

where U(s, 5s) denotes the unitary matrix that diagonahzes the Hamiltonian H(s + 5s) with 
respect to the eigenvector basis of H(s). Thus, U(s) evolves in a multiplicative manner. One 
could think of treating U(s, 5s) by perturbation theory with respect to 5s to derive closed 
evolution equations for the probability distribution of U ka (s) and the corresponding set of 
eigenvalues e a (s). This has, so far, not been undertaken. However, recently Chalker et al. 
f32fl adopted a related picture to investigate the eigenvalue statistics separately by relying 
on simplifying approximations. Here, we only wanted to stress, that the problem of the 
statistics of eigenvectors can be viewed as the problem of a random multiplicative process 
for unitary matrices. 

As to the problem of the statistics of eigenvectors one can easily imagine two extreme 
situations. In the first situation the probability distribution of U is peaked at a single fixed 
matrix. This matrix singles out a certain basis of eigenstates. We choose this matrix as the 
unit matrix, such that the Hamiltonian is diagonal in the initial basis 

H ik = Si5 ik . (69) 

The corresponding eigenstates are localized at certain 'sites', i.e. 

Mk)=5 ak (70) 

The second extreme situation corresponds to an isotropic distribution for the unitary matrix. 
By this we mean that the probability density, V(U), to find a certain unitary matrix U 
within the volume element, d[U(N)], is equal for all elements U G U(N). The volume 
element itself stays invariant under the action of group transformations (invariant measure) 
d[U(N)], i.e. it does not single out any particular element. The corresponding eigenstates 
are also isotropically distributed among all possible eigenstates, no basis is preferred. 

Of course, both extreme situations are not generic ones and we have to see for which 
values of physical scales they may appear. 
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V. SCATTERING MATRIX MODELING 



For mesoscopic conductors with several probes the d.c. conductance coefficients Gkm 



are determined by asymptotic current correlators, Eq. ( |A20| ), and hence for independent 
Fermions at T = by asymptotic Green's functions, Eq. (|A31|) . In a geometry where probes 



are represented as infinite leads attached to the conductor the asymptotic Green's functions 
determine the matrix of transmission amplitudes,^™ > to scatter from channel a in lead k to 
channel (3 in lead m and the Biittiker formula relates them to the conductance, Eq. ( |A36| ) , 



G k ^ m {u = 0, T = 0) = j Tr {uL} (E F ) . (71) 

This formula tells that the dimensionless conductance, especially in a two-probe geometry, 
is the number of effective transmitting modes, 

g = Y, T a, T a = J2\t a p\ 2 , 0<T Q <1 (72) 

a f3 



and it is thus of similar conceptual importance as the Thouless formula, Eq. 

The transmission matrix can be calculated by asymptotic Green's functions. These have 
to be determined under the requirement of attaching the leads to the conductor. For this 
purpose it is advantageous to divide the Hilbert space into two parts by the projection 
operator P onto the Hilbert space of the conductor and its counterpart Q = 1 — P that 
projects onto the leads. Denoting projected operators A by PAP =: A PP ,PAQ =: A P q, 
etc. the general algebraic solution of the problem to calculate the projected Green's function 
reads 

G+ P (E) = (E + - Hpp(E))- 1 (73) 

where the effective (in general energy dependent and non-hermitean) , Hamiltonian Hpp is 
defined via the couplings and the Green's function of the leads Gqq(E) := (E + — Hqq) -1 

Hpp(E) = Hpp + H PQ G + QQ {E)H QP . (74) 

Since Eq. (|73|) is defined for a finite system and the infinite clean leads Green's function 
can be determined analytically for not too complicated geometries the method of effective 
Hamiltonian also helps in doing the calculations numerically. The Green's function evaluated 
at the surface between leads and conductor will then determine the conductance coefficients. 

The Biittiker formula points out that the mesoscopic conductor has very much in common 
with optical wave guides where transmission probabilities between incoming and outgoing 
wave modes are the central physical quantity of interest. It is thus tempting to model 
mesoscopic conductors by scattering matrices in analogy to optical wave guides. 
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A. Quasi-One-Dimensional Conductor 



As a guiding example let us consider a box shaped system (see Fig. 12 with two semi- 
infinite leads attached serving as particle reservoirs (contacts)). The cross section Lf" 1 is 
kept fixed while the the length L is treated variable. Such systems are denoted as quasi- 
one-dimensional (quasi-lD) conductors. The leads are characterized by discrete branches 
of energy dispersion (see Fig. 13) due to the finite extension in transverse direction. For a 
given value of the Fermi energy there exist N c quantum states of incoming and outgoing 
waves (channels). As a rough estimate one may think of N c being the number of lattice 
points one can put on the cross section of the conductor with lattice spacing of half a Fermi 
wavelength, N c « (2L t /\ F ) d - 1 . 

A direct consequence of the Buttiker formula is the integer quantization of dimensionless 
conductance in ideal ballistic conductors where no scattering occurs, or the mean free path l e 
is much larger than L t , L. Then T a = 1 and the conductance displays quantized plateaus as a 
function of the Fermi energy (or particle density). Jumps from one plateau to another occur 
whenever a new channel is occupied, i.e. the Fermi energy crosses a new branch of the leads 
energy dispersion. In contrast to classical physics the conductance of an ideal conductor is 
not infinite but bounded by the finite number of quantum channels. The classical behavior 
can only be recovered for vanishing Fermi wavelength or, equivalently, for infinite number 
of channels. The finite resistance of ideal ballistic conductors can be attributed to the 
contacts where the equilibration of the electrons take place. Therefore, one may think of a 
decomposition of the resistance into a contact resistance, R c = (e/h 2 )l/N c , and an intrinsic 
resistance, R, such that R = R c + R. This yields for the corresponding intrinsic conductance 
9 = Z) a ^a/(1 ~ N' 1 J2 a T a ), a result that was obtained by Langreth and Abrahams |87 



generalizing the one-channel version of Landauer's pioneering works []84| , |85| . However, one 
should keep in mind that the notion of intrinsic conductance is a theoretical construction 
and not a quantity to be measured in a two-probe experiment. Such experiment corresponds 
to the conductance as described by the Buttiker formula. 

In general, the quasi-lD conductor can be described by a 2N C x 2N C scattering matrix S 

(75) 



S 



( rt ,\ 



yt r', 



connecting incoming and outgoing channels left and right of a scatterer, t, r, t' and r' being 
N c x N c matrices of transmission and reflection coefficients for scattering from left to right 
and vice versa, respectively (see Fig. 14). Here t is precisely the matrix of transmission 
amplitudes that appear in the Buttiker formula. Due to probability-flux conservation S is 
unitary, 

S j = S' 1 . (76) 
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To model a long quasi-lD system one can now add several conductors in series. However, 
the total ^-matrix is not multiplicative. Instead of 



(77) 
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= s 











we are seeking a transfer matrix M with the property Af 



1+2 



M 2 Mi, i.e. 



(o^ 




(A 




= M 




v 7 ' 1 







(78) 



A straightforward calculation yields 



M 



\ 



£t) 1 r't'- 1 ^ 



J 



-t'- l r t'- 1 



(79) 



Making use of probability-flux conservation 



(l ^ 
-1 



(80) 



where 1 stands for the iV c x iV c unit matrix we can rewrite the dimensionless two-probe 
conductance as 



g = Tr 



M Aft + (AfAft) 



-i 



(81) 



A modeling of a mesoscopic quasi- ID conductor can now be based directly on the scatter- 
ing matrix or, equivalently, on the transfer matrix. This can be done by fixing the statistical 
properties of the S-matrix corresponding to a small strip of length 5L, denoted as strip 
S'-matrix S(5L) and compose the whole conductor by putting statistically independent strip 
S'-matrices in series. The total S'-matrix follows then from the multiplication of the cor- 
responding strip transfer matrices M(5L). The assumption of statistical independence is 
justified if the strip length 5L is larger than the microscopic potential correlation length ly 
introduced in Sect. |IIB|. Furthermore, this modeling allows for a simple description of the 
mean free path corresponding to the strip S'-matrix. As long as the corresponding reflection 
probabilities | r a p | 2 =: R a p are small compared to 1 the mean free path l e is large compared 
to 5L and can be defined as follows 



^:=iV c - 2 E(^(^)) • 

Le a/3 



(82) 
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The S'-matrix fulfills the requirement of unitarity. To model systems with specific symmetry 
properties one can impose further symmetry constraints on S. The most important case 
is that of time reversal symmetry. This symmetry can be broken by magnetic fields. To 
model systems with magnetic field one faces the problem that the magnetic length (see 
Eq. (|36D) as a chiral length scale has to be incorporated into the model. It is not obvious 
how to do that in the scheme described here. However, weak magnetic fields only break 
time reversal symmetry (the magnetic length being larger than system size) and can be 
simply incorporated in the symmetry properties of S. For example, time reversal operates 
by interchanging incoming and outgoing channels and complex conjugation of amplitudes. 
In case of time inversion symmetry the S'-matrix and transfer matrix M fulfill 



s = s L , m* = >;,..!/>;,. , 



(oi\ 



V 10 / 



(83) 



Internal degrees of freedom, such as spin, can also be incorporated by taken the correspond- 
ing time reversal transformation into account. 

So far, the modeling rests on a microscopic length scale, the mean free path l e , and 
symmetry properties. The statistical problem is thus defined by (i) fixing the distribution 
of the strip S'-matrix S(5L) and (ii) by applying the composition law for the corresponding 
transfer matrices 

M(L + 5L) = M(SL)M(L) . (84) 

This defines a stochastic multiplicative process where the time variable is the system length. 
This resembles Eq. QB8]). (iii) The statistical properties of the increment 

5M(L,5L) = M(L + 5L)-M(L) (85) 

are known by construction. Following the steps (i) - (iii) it is possible to construct a stochas- 
tic differential equation (Fokker- Planck equation) for the distribution function V(M; L) (see 
e.g. |0]). Still, a general solution of this equation is presently not available and we will not 
discuss it in detail. A related approach, based on simplifying assumptions, will be discussed 



in Sect. VIII B 



The matrices ttf and MM^ occurring in the conductance formulas, Eqs. ([71], |8T[), are 
hermitean and, thus, can be diagonalized. The positive eigenvalues are denoted as < % < 1 
(for ftt) an d as < < oo (for MM f ). The eigenvalues of MM t appear m inverse pairs 
and we can restrict to those with Vi > 0. The conductance reads in the corresponding 
eigenvalue representation 

» = ?* = ?I+JiM (86) 
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Here the eigenvalues % have to be distinguished from the total probability to be scattered 
from a channel a into forward direction, T a , that occurs in Eq. flT2"|). A related polar 
parameterization (cf. ||121|| ) of the transfer matrix M in terms of radial coordinates < Aj < 
oo and four unitary iV c x N c matrices u l is also used frequently 



M 



u 1 
u 3 



Vx v / TTA 



W 
u A 



(87) 



Here A = diag(Ai, . . . , AjvJ and 



It worth noticing that in these S'-matrix models the conductance itself appears as a linear 
statistics | of eigenvalues of a certain random matrix. However, one should keep in mind 
that the statistical properties of these eigenvalues are, in general, not independent of those 
of the associated eigenvectors. 



B. Localization 

Before we introduce generalizations of the S'-matrix modeling we will show that the 
present one is well suited to discuss the problem of localization in quasi- ID conductors. 

An instructive example is the ID conductor for which N c — 1. The composition law 
Eq. flgip tells that the composition law for the transmission reads (see Fig. 15) 

tn = • (89) 

1 - r[r 2 

Exploiting unitarity of S yields for the transmission probability T — |t| 2 = l — R, R = \r\ 2 = 
r\2 



\r 



Ti 2 = =^L= (90) 

1 - 2cos0 v / ^i^ + J Ri^2 



where is the sum of phases of r[ and Based on Eq. (|90|) one can derive an evolution 
equation for the probability distribution of T with increasing length L. We will postpone 
this to Sect. [VII C| . Here we will take advantage of the fact that an average of lnT 12 over a 



uniform distribution of phases <fi yields the simple expression 



3 By linear statistics one denotes quantities that are sums of a specific function / of the random 
parameters, J2if(^i)- 
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<lnT 12 )^ = lnTx + lnT 2 . (91) 

Thus, the phase-average of the logarithm of the transmission turns out to be additive. 
Iterating this procedure tells that the logarithm of the transmission will be distributed in 
a Gaussian way according to the central limit theorem for independent additive random 
numbers. This motivates to call T t := exp (lnT) the typical transmission, and to write down 
a quantum series composition law for conductance g = T in ID conductors 

g t {L 1 + L 2 )=g t {L 1 )g t {L 2 ) (92) 

This composition law is in striking contrast to the macroscopic composition law for incoher- 
ent conductors, G _1 (Li + L 2 ) = G~ l (L\) + G _1 (L 2 ). Equation (|9~2^) immediately tells that 
the conductance will exponentially decrease with increasing system length L, which means 
localization. To determine the localization length £ we calculate 

d 9t(L) fT ^ 9t(SL) - 1 

-ir = 9 ^ L) 6L ' 5L ^°- (93) 

Thus, g t (L) = g t (L ) exp(-2(L - L )/£) with £/2 = 5L/{1 - T t (8L)), and by Eq. © 

£ = 2Z e . (94) 

This result demonstrates that in one-dimensional conductors all states are localized. Fur- 
thermore, the localization already occurs on the scale of the mean free path. This is a crucial 
mesoscopic effect as it is a consequence of the multiplicative quantum series composition, 
Eq. ©. 

For arbitrary channel numbers N c the above ID result does not apply. However, one can 
make an educated guess what might happen to the localization length, if localization occurs 
at all. 

To this end, we notice that the conductance in ID for length L much smaller than l e 
is given by g(L) = 1 — L/l e in accordance with our definition of l e . Since the result g = 1 
corresponds to ideal transmission of the wave guide an Ohmic series decomposition into 
contact resistance 1 and intrinsic resistance, L/l e , is meaningful. Assuming that for N c > 1 
and L £(iYc) Ohm's parallel composition law is applicable yields the classical conductance 

9d(L) = ^ (95) 

and the natural limit for this law to hold, g(L rs £) w 1, yields an estimate of £(N C ) 

£(AQ w iV c / e . (96) 

We are now going to discuss rigid results that confirm the localization effect for arbitrary 
iV c . In Sect. |V111B| we will see that our estimate of the quasi-lD localization length is in 



accordance with known results. 
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A crucial point in the ID calculation was the multiplicative composition law for typical 
conductance. The transfer matrix obeys also a multiplicative composition law and the 
question arises what can be said about the statistical properties of a large product of random 
transfer matrices. Recall that the product of independently distributed positive random 
numbers X(i), 



N 



X(N) = l[X(i) = exp 



i=i 



' N 
i=l 



(97) 



gives rise to a central limit theorem for its logarithm: For large N, the distribution of 
In X(N) is well approximated by a Gaussian with mean value N (\nX) and variance 



a 



N 



N 



(lnX) 2 ) - (\nXy 



(98) 



For the original variable X(N) the distribution is called log-normal distribution and reads 



V(X,N)dX 



2txo~ n 



exp 



2a% 



(InX -N {lnX)f 



dlnX 



(99) 



Furthermore, one can conclude that for almost all realizations the geometric mean asymp- 
totically coincides with the typical value 



lim {X{N)) 1/N = exp ((lnX)) =: X t . 



(100) 



There exist some extensions of the law of large numbers to the case of products of 
independent random transfer matrices (see ||37|| ). A version of the theorem of Oseledec | 103| 
guarantees the existence (with probability 1) of eigenvalues of the diagonalizable limiting 



matrix 126,127 



M := lim (M^(N)M(N))^ 



(N = L/Lq, Lq some fixed reference length) of the form 



M 



"7i 



)■ 



(101) 



(102) 



Motivated by the analogy between transfer matrices and time evolution matrices considered 
in chaotic dynamics the quantities r y 1 < . . . < j N are called Lyapunov exponents. 

The Lyapunov exponents 7, are asymptotic non-random numbers and are the so- 
called self averaging limit of the random numbers Ui(N), the exponents of eigenvalues of 
M\N)M(N) 1 divided by 2N 



li 



lim 



N->oo 2N 

A system with channel number iV c has thus N c characteristic length scales 



(103) 
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6 = - = lim ( 104 ) 



According to the conductance formula Eq. ( jS6|) the inverse of the smallest Lyapunov 
exponent yields the relevant localization length 

€ = - • (105) 

7i 

A theorem by Furstenberg states that the largest Lyapunov exponent is always 
positive. This proves localization in strictly ID, but is meaningless for iV c > 1. Virster ||127| 



has later proofed that j t is strictly positive, as long as N c is finite. Thus, in the quasi-lD 
limit L — > oo the localization length £ is finite. 

We postpone the discussion of the statistics of Ui(L) to Sect. [VIII B| and proceed by 
introducing more general S-matrix models. 



C. Networks Of Scatterers 

So far, the modeling of mesoscopic systems by scattering matrices was restricted to the 
quasi- ID case. More general c?-dimensional mesoscopic systems can be described by a model 
introduced by Shapiro [ |114|J44]| . Instead of taking a strip S-matrix one can consider a regular 



network (lattice spacing a) consisting of N = volume/a d sites and Nj, bonds (see Figs. 16, 
17). Each site represents a scatterer and is described by a unitary 2d x 2d scattering matrix 
S. The bonds represent free propagation between the blocks. Each bond carries two waves 
propagating in opposite directions. Each S'-matrix transforms 2c? incoming amplitudes into 
2c? outgoing amplitudes. To each bond amplitude a random phase is added which simulates 
an irregular lattice of scatterers which are distributed between distance a and a + Xf where 
Xp is the wavelength of scattered waves. 

By specifying appropriate boundary conditions one can describe multi-probe conductors 
within this network model. Consequently, physical quantities like conductance coefficients 
and even channel dependent transmission amplitudes can be studied with less effort than in 
systems described by a Hamiltonian. 

Interesting is the fact that the S'-matrix modeling does also allow to determine energy 
eigenstates and eigenvalues. This comes along with the following interpretation: mapping 
incoming to outgoing states by scattering matrices usually refers to a process consuming 
infinite time in order to guarantee energy conservation. However, if one allows for a certain 
uncertainty in the energy, a finite time is enough to follow the scattering process. Therefore, 
one can take the point of view that the scattering matrices at each site describe a unitary time 
evolution of states in one unit of time. By this interpretation the finite energy uncertainty 
is treated as negligible. Accordingly we define the unitary matrix of time evolution in a unit 
of time by [H,H,^1 
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J2U ik ^(k;t)=^(i;t + l). (106) 
k 

Here {ip(k;t)} forms, for each time t, a set of 2Nt> complex bond amplitudes. Normalizing 
them, Y^i IVKO 1 2 = I; makes them a proper choice for local wavefunctions. The unitary 
2Nt> x 2Ni, matrix U is uniquely determined by the network of scattering matrices since it 
summarizes nothing but the scattering conditions at the sites. Energy eigenstates ip are 
stationary solutions of the unitary operator U 

Uip^ip (107) 

and thus correspond to the eigenstates of U with eigenvalue 1. The energy of the underlying 
electron system enters here only parametrically through the scattering amplitudes. For 
simplicity, let us assume that the scattering strengths are equal at all nodes (abbreviated by 
one symbol s) and randomness enters only through the bond phases. Then U will have an 
eigenstate to eigenvalue 1 only for a discrete set of scattering strengths, which correspond 
to a discrete set of energy eigenvalues of the underlying electron system. Instead of looking 
for the discrete set of energies one can also fix the energy and ask for the more general 
eigenvalue problem 

U{E)^ a {E) = e iME) ^ a (E) • (108) 

ip a (E) are, for fixed E, a set of 2Nb eigenstates at energy E of slightly modified disorder 
configurations, each of the modifications being an overall shift in random phases on the 
links. The eigenphases 4> a (E) can be considered as an energy excitation spectrum at energy 
E. Thus, the network models allow to investigate energy eigenvector statistics and, via the 
phases <fi a (E), energy eigenvalue statistics. 

The network models have a further advantage as compared to Hamiltonian models con- 
cerning numerical calculations. In case that only a finite resolution of energy scales is 
required one does not need to fully diagonalize the unitary operator. Instead, one iterates 
the time evolution step by step until a certain time scale is reached, the inverse of which 
yields the required energy resolution. 

Finally, we like to point out that within the network models there are possibilities not only 
to include global symmetries but also characteristic length scales. For example by allowing 
for internal degrees of freedom (e.g. spin) the number of amplitudes on each bond can be 
increased and scattering strengths can be introduced that take the change of internal degrees 
of freedom into account. The inclusion of strong magnetic fields with their characteristic 
chiral character is also possible. No commensurability features like in the Anderson model 
with Peierls substitution do occur. Especially, one model of this category has become widely 
known as the Chalker-Coddington model in 2D which only allows for scattering to the 
left and right (see Figs. 18, 19). 
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VI. STATISTICS IN IDEALIZED SITUATIONS 



In this section we consider idealized mesoscopic systems and ask for the statistics of 
energy eigenvalues and eigenvectors. The idealization is due to the assumption that the 
statistical properties of both can be studied separately. We distinguish between the local- 
ized situation imitated by a diagonal Hamiltonian and the delocalized situation imitated 
by a standard random hermitean matrix having independent Gaussian distributed random 
entries. 

We will review the energy level statistics in these systems and the statistics of the local 
box-probability P(L b ) for an electron that is in an eigenstate ^(r) to be found in a certain 
box of volume Lf, 

P(L b ) := f d d r | ^(r) | 2 . (109) 



box 



The complete information about the (energy) level statistics of a N x N Hamiltonian is 
contained in the joint probability density 

V(e 1 ,...,e N ) . (110) 

To simplify notations we will assume throughout that any reordering of levels in V does not 
change the probability. The simplest statistical quantity to be extracted is the average level 
density 

u(E) := (p(E)} = N- 1 (^2 S(E - = j de 2 . . . de N V (E, e 2 , . . . , e N ) (111) 

defining the average level spacing A(E) = (Nu(E))~ 1 . The two-level correlation function is 

R(s)= A 2 (^5^-6^5(^-6^- I, s := . (112) 

For simplicity we assumed that it is translational invariant within a band of states with 
A being independent of energy. By writing the energy separation in units of average level 
spacing (s) one focuses on those properties that are independent of the chosen energy units. 
Quantities which contain information about n-level correlation functions are the n-level 
spacing distribution Q(n, s) which yields the probability to find exactly n levels in an energy 
interval of width s (the center of which can be arbitrary). The probability to find exactly no 
level inside the interval defines the so-called level spacing distribution P(s) = d 2 Q(0, s)/ds 2 
which yields the probability density to find a level spacing s with respect to ds. The average 
number (n(s)) of levels in an interval s and its variance £ 2 := (n 2 (s)} — (n(s)} 2 are also of 
interest. The derivative of £ 2 with respect to (n) defines (for large (n)) a level compressibility 

dT? . . 

X ■= . \im -jj-r ■ 113 
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A. Localized Phase 



We consider iV orthogonal energy states ip a , a = 1, 2, . . . , N each of which is localized 
to one of iV sites of a regular lattice (lattice spacing a) in an arbitrary dimension d, leading 
to a vanishing localization length. The corresponding energy eigenvalues are independent 
random variables taken from a distribution 



e 2 



P(e) = -«=-exp(- s3 . (114) 



V2TTE0 V 2^0, 

The corresponding random Hamiltonian is diagonal in site representation 

H ik = £i$ik (115) 

with eigenstates 

ip a (i) = 5 ai for e a = (116) 
This situation was already addressed in Sect. |IV C| as one of two extreme cases for the 



eigenvector statistics. As compared to the Anderson tight-binding Hamiltonian, Eq. ([£]]), the 
Hamiltonian in Eq. (|115|) ignores the kinetic energy completely. The energy level statistics 
is contained in the joint probability distribution 



N 



V(e 1 ,...,e N )= l[p(e a ). (117) 



Of=l 



Consequently, the average level density is 

u(E)=p(E) (118) 



and the average level spacing within the band (\E\ Eq) is A = l/(Np ), p := (v27T£o) _1 - 
The two-level correlation function R(s) reflects the fact that there is only auto-correlation 
between levels (up to 1/N corrections). 

R(s) = 5(s) - N- 1 (119) 

The level spacing distribution follows the Poisson law 

P(s) = exp(-s) (120) 

(that is why the ensemble is called Poisson ensemble) and the characteristic relation between 
number variance and average number of levels can be obtained 

S 2 = (n) (121) 

indicating that the levels are compressible with compressibility 1. 
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The probability to find an electron in a box of size a in the eigenstate ijj a is 

P a (a;i) = \^ a (i)\ 2 = 6 ai . (122) 

Increasing box sizes to L b > a does not change the result as long as the boxes are chosen 
to be non-overlapping and centered at the original point. Independently of the box size the 
probability is either zero, if the state is not located in the box or 1 in the opposite case. 
For a given energy E we consider the microcanonical average of P a 

P(E,L b ;i) ^W-e«)P«(L b -,i) 

Y,pd{E - Bp) 

where a smearing-out of 5-functions is understood. The smearing width T(E) has to be 
adjusted such that isolated 5-peaks are broadened, but not over many neighboring levels. 
A convenient choice is to take T(E) of the order of the actual level spacing at E in a 
given realization. The denominator in Eq. (|123|) is then 1/T(E). Since the eigenvalues 
are independently distributed, the centers of the corresponding eigenstates are randomly 
distributed over the lattice points. Therefore, in almost any case one will find zero box- 
probability at a certain site i, unless the energy E coincides (within the scale of the smearing 
r(i£)) with eigenvalue e$, in which case the box-probability is 1. The average box-probability 
is finite and given by (P(E)) Lb = (l/N)(L b /a) d = {P) w 

Thus, the probability distribution V(P; L b ) of the box-probability is obtained as 

V(P; L b )dP = \8{P) (l - (P) L ) + 8{P - 1) (P) L ] dP . (124) 



The average value is not characteristic for this distribution, since the main weight is given 
to P = 0. However, e.g. the geometric mean 

P t (L b ) :=exp(lnP) Lb = (125) 

reflects a typical box-probability. Alternatively, one may consider the median of the distri- 
bution (that value up to which half of the total weight is used up) which vanishes also and 
could be termed typical as well. 

The LDOS p(E,i) = ^2 a 5{E — e a )\ip a (i)\ 2 with a smearing-out of 5-functions on the 
scale of T(E) yields for the LDOS in a box, p(E, L b ; i) 

p(E,L b ;t) :=T{E)- l P{E,L b -i). (126) 

However, we must mention that this notion of LDOS in a box differs from that used by other 
authors (see e.g. [0). If one considers the LDOS with a parametric Lorentzian smearing-out 
of 5-f unctions, 

p(E;r)=Y— \ M«(r)| 2 , (127) 

^ V ; V27r(£-£ a ) 2 + r 2 /4^ V ; ' V ; 
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where V/A can vary from to oo one can investigate parametrically the interplay of level 
and eigenvector statistics. 

To simplify the discussion we will be often more restrictive and call 

p{E,L b ;i) := A^PiE.U-i) (128) 

the reduced LDOS in a box. This definition concentrates on the eigenvector statistics. Con- 
sequently, the distribution function reads 

V(p; L b )dp = [8{p) (l - A (p) Lb ) + 5(p - A' 1 ) A (p> J dp . (129) 

In summary, ideal localized systems (with vanishing localization length £) are character- 
ized by a compressible spectrum of uncorrelated energies and a LDOS in a box that typically 
vanishes, although the average DOS is finite. 



B. Delocalized Phase 

For the idealized delocalized phase we make the assumption that eigenstates are dis- 
tributed isotropically within the space of all possible eigenstates for a N x N Hamiltonian 
matrix, Hi k . In other words the diagonalizing unitary matrices U ka , Eq. (|65|), are uniformly 
distributed with respect to the invariant measure d[U(N)]. This situation was already ad- 
dressed in Sect. [IV U| as the second extreme case for the eigenvector statistics. 



Therefore, we assume the distribution of the Hamiltonian matrix to be unitarily invariant. 
One of the simplest possible choices is the so-called Gaussian unitary ensemble (GUE) 
introduced by Wigner ||134j| when studying the level statistics of complex nuclei, 



V ({Re H ik , Im H lk }) d [H] = C N exp ("^Tr H 2 ^j d [H] . (130) 

Here CV is a normalization constant, Sq some arbitrary energy scale and the volume element 
is defined in terms of the independent matrix elements of H as 

N N 

d[H] = l[dH ll l[d(Re H lk )d{lm H lk ) . (131) 

1 i<k 

Since Tr H 2 = 2 J2 (R- e H ik ) 2 + (Im H ik ) 2 + J2iHfi the GUE describes a random matrix 

i<k L J 

with all its elements uncorrelated. Each of the elements vanishes on average and its absolute 
value fluctuates, (\H ik \ 2 ) = Sq/N. If we think of the model in site-representation and 
compare it to the Anderson model we see that by Eq. ( |130|) hopping terms are dominant 
and hops from one site to another are equally probable independent of the spatial distance. 
This means that we can associate a vanishing time scale tr> to travel through the system. 
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Introducing eigenvalues e a and eigenvectors 4> a (k) = Uk a one can transform the proba- 
bility V({Re H^, Im Hik})d [H] to these variables at the expense of introducing a Jacobian 
between the set of variables. Fortunately, the Jacobian can be calculated, the integration 
over d\U(N)] being a trivial normalization (since Tr H 2 is unitarily invariant), and one ob- 
tains for the joint probability of eigenvalues 

/ N N \ 

V(e ll ...,e N ) = C N l[(e a -e ) 2 exp(-J2—^e 2 a ) . (132) 

a<(3 \ a ^0 / 

The factors in front of the exponential are due to the Jacobian and ensure that the prob- 
ability to find two levels close to each other vanishes; a phenomenon which is denoted as 

level repulsion. The factors can be rewritten as exp 2 hi \e a — £p\) such that the joint 

probability describes a classical Gibbs ensemble, 

V ({e a }) = C N exp [-/3H ({£„})] , (133) 

of a gas of particles with coordinates e a and a Hamiltonian 

W ({*«}) = \H U(e a , e p ) + £ V{e a ) (134) 

that contains a logarithmic two-body interaction U(x,y) = — In |x — y\ and a one-body 
(confining) potential V(x) = Nx 2 /(2(3£ 2 ). The inverse temperature is (3 = 2 for the GUE. 
A related ensemble of real symmetric matrices, reflecting time inversion symmetry, denoted 
as Gaussian orthogonal ensemble (GOE), gives rise to the same Gibbs ensemble, the only 
change being /3 — 1 (/3 = 4 corresponds to spin systems with time reversal symmetry, for 
more details and review see The interpretation of the joint probability distribution of 



levels as a Gibbs-ensemble has led to an important mean-field approach to level statistics 
43] which we will use later in the context of quasi-lD conductance (cf. Sect. |VTI FBI). In the 



standard random matrix of GUE and GOE a number of results are well known: 
1. The average DOS is (for N ^> 1) given by Wigner's semi-circle law, 

1 l-f^-V. (135) 



with average level spacing (E ~ 0) A = tt£ q /N. 

2. The two-level correlation function shows level-repulsion for s<l and decays ~ — 1/s 2 
for s > 1. In the GUE it reads (s ^ 0, iV -> oo) 

, , /sin7rs\ 2 
R(s) = -( . 136 

V 7TS J 
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3. 



The level spacing distribution is very well approximated by Wigner's surmise 

P(s) = Ap s p exp {-Bp s 2 ) (137) 



where e.g. A 2 = tt/2,B 2 = <i/ir. The characteristic small s behavior ~ s 13 reflects the 
power of level-repulsion characteristic of the symmetry of the Hamiltonian. 

4. The level number variance £ 2 depends only logarithmically on the average level num- 
ber, 

£ 2 oc In (n) , (138) 

which yields a vanishing level compressibility (~ (n)" 1 — > 0). 

Eigenvectors ip a {k) = Uk a are isotropically distributed uncorrelated Gaussian variables 
up to l/N corrections. More precisely, in the GUE the eigenvector components are complex 
with uncorrelated Gaussian distributed real and imaginary parts. The average of each is 
zero while the variance is, due to normalization, given as 

((Re Mk)f) = ((Im ip a (k)) 2 ^ = ± . (139) 

For the GOE the components are real with vanishing average and variance {\^ a (k)\ 2 ) = N" 1 . 

This fixes the distribution of the probability P(E; a) to find an electron in a box of size 
a (a lattice spacing). In the GUE one finds the Rayleigh distribution 

V(P; a) = N exp (-PN) , (140) 

and in the GOE the Porter-Thomas distribution 

V{P]a) = \2^p) ex P(" PiV / 2 )- ( 141 ) 

Both distributions follow immediately once an intensity of waves ip is considered which are 
large sums of either complex or real numbers. This leads, by the law of large numbers, to 
Gaussian distributions for ip. Both distributions can be generalized to arbitrary box size 
L b by replacing iV" 1 with (P) L = A^ _1 (L; ) /a)' i . Characteristic of these distributions is the 
exponential tail and the fact that no parameter besides (P) Lb occurs. Still, the average value 
is a typical value for the distribution; e.g. the median differs from the average value by only 
a factor of In 2. 

The statistics for the reduced LDOS in a box, p(E; L b ) = A~ 1 P(E; L b ), is then straight- 
forward. 

In summary, ideal delocalized systems (with vanishing time scale to) are characterized 
by an incompressible spectrum of correlated energies and a LDOS in a box the distribution 
of which develops an exponential tail. Still, the average value is typical. 
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VII. CONCEPT OF SCALING 



The level and eigenvector statistics of the idealized systems show a remarkable feature: 
No physical scale besides the kinematic scales, level spacing A, total volume Na d and box 
volume Lf enters the statistical quantities. The parameters of the models, Eq (for the 
localized phase), and So only define the energy unit for A. No transport scale enters the 
models. This is consistent with the following interpretation: In the localized phase the 
only relevant scale is the localization length, which was treated as zero. In the idealized 
delocalized phase the only relevant scale is the time to cross the system, or equivalently 
the corresponding Thouless energy E^h = h/tp which was treated as infinite. We may 
thus expect that the statistics described in the foregoing section correspond to asymptotic 
situations for which the dimensionless conductance is either zero (localized phase) or infinite 
(delocalized) phase. 

The question then arises how the statistical properties do change for finite conductance. 
Of particular importance are the statistical properties when approaching the LD transition. 
Which physical scales determine the statistical properties? A related question concerns the 
change of the statistical properties with increasing system size. 



A. Relevant Scales And /3-Function 

Looking at the conductance formulas by Thouless and Biittiker, Eqs. ( |20| , |72|) , it is tempt- 
ing to assume that the conductance itself is the (perhaps only) relevant quantity to distin- 
guish localized and delocalized phases and to characterize the LD transition. 

1. The localized phase with 

9 « 1 (142) 

corresponds to only closed transport modes (channels) (Biittiker formula with T a <C 1) 
and the existence of a finite localization length £ <C L, g ~ exp(— 2L/^). 

2. The delocalized phase with 

9 » 1 (143) 

corresponds to many open channels (Biittiker formula with N c m g modes of T a m 1) 
and a Thouless energy much larger than the level spacing, S> A. 

3. The LD transition is expected to occur when scales match, i.e. 

E Th ~ 0(A) , g~ £>(!) . (144) 
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According to the Biittiker formula only a number of (9(1) modes are open and the local- 
ization length diverges. This divergence turns out to be algebraically on approaching 
the transition point, 



— v 



(145) 



where r measures the deviation of a tuning parameter (e.g. Fermi energy) from its 
critical value at the LD transition. The exponent v is called critical exponent of the 
localization length. 

By this reasoning a relevant length scale £ c can be defined via the conductance: £ c is 
that fictitious system size for which 



In the localized phase £ c can therefore be identified with the localization length, while in the 
delocalized phase £rh(£c) ~ C(A(£ C )). 

As our considerations call for a description of the LD transition as a critical phenomenon 
we briefly review some aspects of the latter. 

A thermodynamic system is described in terms of the coordinates of its equilibrium states 
manifold M. (e.g. volume V, particle number N and temperature T). Any thermodynamic 
quantity of interest can be calculated from a thermodynamic potential (e.g. the free energy 
!F(y,T,N)) which is an analytic function of these coordinates (e.g. the specific heat cy = 
d 2 J r /dT 2 ). Generally, one can change the parameterization of M (e.g. volume V — > pressure 
p) to describe the same equilibrium state - at least locally. In the thermodynamic limit 
where N and V tend to infinity simultaneously leaving the density n = N/V constant, 
it may happen that the thermodynamic potential is no longer an analytic function of its 
coordinates. In such a situation the state of the system is not uniquely described by only 
one set of coordinates manifesting in phases. When crossing the coexistence curve the 
thermodynamic potential jumps. Such a process is called a first-order phase transition. 
The jump is quantitatively described by the value of the latent heat. Assume that the 
coexistence curve ends in some point of the (p, T) diagram which is called the critical point. 
Here the latent heat vanishes, the thermodynamic potential is continuous and the difference 
of densities between the phases vanishes. Yet, at the critical point, the thermodynamic 
potential is not analytic since its second derivatives are discontinuous. It is then helpful 
to look for a quantity which serves as an order parameter of the transition. It should be 
finite in one phase and vanishes in the other phase. A more restrictive requirement is that 
an order parameter should correspond to a local quantity <p(r) the average value of which 
m = (y?(r)) vanishes with some power law at the critical point (e.g. the difference of the 
corresponding densities in the transition), 



<7&) ~ O (1) . 



(146) 
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m = , T <T C 

mocr?, T>T C . (147) 



The exponent ft is called the critical exponent of the order parameter and r = \{T — T c )/T c \. 
It is a necessary condition for such power law behavior that the phase transition is not of first 
order since a power law is the paradigm of a physical law with no scales involved. In a first 
order transition, however, the jump in the thermodynamic potential sets a physical energy 
scale. A defining characteristic of a critical phenomenon is that there is also no relevant 
length scale at the critical point. Thus, the correlation length £ c of the order parameter, 
m = ((p(r)), defined by the statistical properties of the local order parameter field 

X {r) := (<p(0)<p(t)) oc exp (-r/f c ) , r := |r| , (148) 

has to diverge at T c with a power law with some critical exponent u, 

£ c oc r~ v . (149) 

The crucial assumption in any scaling approach to critical phenomena is that the critical 
exponents ft, v have their origin in the divergence of a single relevant length scale, the 
correlation length £ c . At the critical point where £ c diverges the correlation function obeys 
a power law 

x(r)ocr~^. (150) 

Using £ c as a cutoff when calculating the global susceptibility x — j d d r x{ r ) we find a scaling 
relation 

7)=—. (151) 

Let us now return to finite systems of linear size L and assume that both the distance 
r and the finite system size L are within the regime without length scales and that the 
function C(r,L) := ({p(0)ip(r)) L obeys a power law with respect to both lengths 

C{r,L) (xr~ v L' y . (152) 

The new exponent y describes the system length dependence and r\ ^ fj describes the distance 
dependence. To consider this case turns out to be relevant in the LD transition problem as 
will be discussed in more detail in Sect. |LX A|. A similar reasoning as above yields now 



y + r\ — i] — — . (153) 

Notice that the quantity r\ is not determined by ft and v alone but requires the additional 
knowledge of the unusual exponent y. 
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Far away from the critical regime there may still exist a scaling variable A(L) which 
replaces the role of T in the following sense. The state of the system is described for a finite 
size L already by the quantity A(L) and the flow with increasing system size L is determined 
by a single function of this variable, the so called /3-function 

P { lnA):= d J^l (154) 
which should be a function of In A alone. Provided the /3-function is smooth, the flow is 



determined by the solution of the differential equation, Eq. (|154| ), frequently called a renor- 
malization group equation. The regime where /3(lnA) can be linearized around a fixed point 
A* of the flow (/?(lnA*) = 0), 

(3 (In A) = 0' (In A - In A*) , (155) 

with f3' being the slope of the /3-function at A*, is called the critical regime. It is in this 
regime where power law scaling can be expected. If we start with a system of size L = L , 
A = A chosen close to A*, and turn on the renormalization flow until the system reaches 
a size £ c determined by the range of validity of Eq. ( |155| ) (which is here by definition the 
correlation length) we get 



/lnA-lnA* 



L J V In A - In A* 



(156) 



Assuming that the critical regime is narrow, we can expand A(£ c ) and Ao around A* and 
get for the correlation length 

/ A — A* \ 1//3 ' 

e e =Mi — -) . (157) 



jYo-A*, 

The width of the critical regime AA = A — A* is determined by r such that 

£ c oc (AA)~ V oc t~ v . (158) 
Thus, the exponent v is given by the inverse slope of the /^-function at the fixed point 

v = l/j3'{h*). (159) 
The behavior of the scaling variable close to the transition is therefore of the form 

A(L; r) = A* + AtL 1 '" + O (r 2 ) . (160) 



The coefficient A can be brought to the form AtL 1 ^ = A (L/£ c ) 1 ' 1/ where A is of order A*. 
Note, that, in general, the coefficient A can differ on both sides of the transition, while the 
critical exponent v is the same. 
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It is crucial that any other scaling variable A which is a smooth function of the previous 
one leads to the same value of v, although the overall form of its /3-function can differ. This 
reflects the fact that critical exponents are universal numbers characterizing a whole class 
of systems. Such systems may differ by microscopic details, however the critical behavior is 
characterized by universal critical exponents of a few relevant quantities. 

In Sect. |1V B| we have learned that knowing the Green's function is, in principle, enough 
to describe the transport properties including the LD transition. A natural guess for the 
order parameter field is the LDOS since it typically vanishes in the localized phase and is 
finite in the delocalized phase. Consequently, the average DOS (which is determined by 
the one-particle Green's function) appears as the natural (formal) order parameter. Since 
the average DOS does not show the LD transition but is a smooth function of the energy 
(which, in the problem of the LD transition, is the analog of temperature in equilibrium 
phase transitions) the critical exponent (3 of this formal order parameter vanishes 

p = 0. (161) 

The LD transition is not manifested in the one-particle Green's function, but in transport 
related quantities such as the two-particle Green's function. We like to mention that this 
fact rules out the possibility of a straightforward mean-field theory for the LD transition. 

Furthermore, as the scaling relation Eq. ( |151| ) tells us the correlation exponent of the 
local susceptibility (a two-particle Green's function) vanishes, 

z = 0. (162) 

This result seems to rule out a power law behavior with finite power for the two particle 
Green's function (say the density correlator) at criticality. In spite of that a finite power has 
been observed []27| for the density correlator at a LD transition. However, it turned out that 
it corresponds to a regime where the correlation function shows a power law (with unusual 
exponent y) also with respect to the system size L. 

A further unusual feature of the LD transition, which at first glance seems to be inde- 
pendent of the previous ones, was first pointed out by Al'tshuler et al. |J when reexamining 
the phenomenological scaling theory of Abrahams et al. fl|]. 

This theory picks up the ideas that we introduced in the beginning of this section. The 
conductance g(L) of a cube with linear dimension L is treated as a scaling variable, i.e. g 
obeys a differential equation in terms of a /^-function which, for a certain universality class, 
is a unique function of Ing. The fact that one scaling variable is enough to describe the 
transition is referred to as one-parameter scaling theory. 

That such treatment makes sense can be seen when calculating the /3-function for extreme 
situations. In a good metal with a conductivity that is a material constant, Ohm's law, 
Eq. ([17D, yields 



P(\ng)=d-2 (163) 
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while in the exponential localized phase one obtains 



(3(lng) = lng + A d (164) 

where Ad is a constant for each dimension d. Assuming (3(lng) to be a smooth function 
one may interpolate between the asymptotic branches (see Fig. 20). From this picture one 
concludes: In ID all states are localized, in 3D a generic LD transition takes place and for 
2D the situation depends on the sign of weak-localization effects for large g. If these tend 
to localize, the /3-function stays negative and all states will localize in the large L limit. 
For weak anti-localization (found in 2D systems with spin-orbit interaction p3| , |l8| ) a LD 
transition should occur. It has been observed by numerical model calculations (see e.g. [|5T|j). 
In an experiment by Kravchenko et al. f8(| on a strongly interacting 2D electron system a 
LD transition was observed. At present it is not clear if it can be due to strong spin-orbit 
scattering. The quantum Hall effect, as mentioned in the introduction, is exceptional in 
this picture. This is due to the absence of a metallic phase. The corresponding /^-function 
terminates from below at the critical point. 

Although the one-parameter scaling theory by Abrahams et al. [p] is simple, predictive 
and qualitatively correct, it has an obvious shortcoming as was pointed out by Al'tshuler et 
al. || . The conductance of a mesoscopic system strongly depends on the individual system 
properties (e.g. on the given disorder potential). As a result, an ensemble of different sys- 
tems exhibits conductance fluctuations which can become so strong that the corresponding 
distribution cannot be characterized by the mean value alone. The latter as well as higher 
moments, can be drastically influenced by the far tails of the distribution. It can thus 
happen that the mean value of the conductance (g) can not serve as a scaling variable. 

Instead, an appropriate scaling approach to the LD transition has to be set up in terms of 
the complete distribution function of the conductance. However, as pointed out by Shapiro 
[ 116| it may be possible to apply the renormalization group approach to certain parameters 



occurring in that distribution function. If an appropriate parameter of the distribution exists 
which is less influenced by the far tails and which can serve to define a typical conductance, 
g t , it may be possible to find a flow equation in terms of a /3-function for this typical 
conductance. In this sense the theory by Abrahams et al. can still be complete. For this 
scenario to be realistic, at least in the vicinity of the transition point, only one length scale 
£ c should diverge at criticality. Then the LD transition obeys one-parameter scaling and the 
exponent v of the correlation length £ c is given by the /3-function of g t . 

It is worth mentioning that the conductivity, defined by the conductance as cx(L) = 
G(L)L 2 ~ d , is often referred to as an order parameter of the transition. It vanishes at the 
transition point in d > 2 and a(L — > oo, r) oc r s defines an critical exponent s. However, 
this notion of conductivity gives no further information than that already contained in the 
conductance. Only for d > 2 the conductivity does vanish at the transition point and the 
exponent s is trivially related to v by s = (d— 2)u since G(£) — G* oc t^}I v . This relation is 
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sometimes referred to as a scaling relation between critical exponents s and v suggesting that 
this replaces the scaling relation Eq. ( |151| ) in common critical phenomena. To the author this 
interpretation is misleading. Firstly, the conductivity is defined via the conductance, which 
is the scaling variable of the problem. For any critical phenomenon with a scaling variable 
A(L) one can define similar order parameters by simply multiplying the scaling variable 
by a negative power of L. It is obvious that this does not give further insight. Secondly, 
for LD transitions occurring in d = 2 the conductivity stays finite at the transition, which 
invalidates the notion of an order parameter for a. Thirdly, the conductivity is not related 
to a local field of the problem. 

In Sect.|X| we will see that the choice of an appropriate order parameter and the occur- 
rence of broad distributions of physical quantities are intimately interrelated. 



B. Scaling Of Typical Values 

In this section we will exploit our experience with the notion of typical conductance 
in ID and make an extension to higher dimensions following a work by Shapiro [TTB ]. The 



extension relies on uncontrolled approximations and may thus be only of pedagogical interest. 
However, the results are surprisingly far-reaching and consistent with commonly accepted 
results about LD transitions. 

Recall that the typical transmission probability in ID followed a simple multiplicative 
composition law which under a scale transformation L — > bL leads to the scaling law in ID 

T t (bL) = {T t {L)) b . (165) 

We would like to generalize this law to higher dimensions. The multiplicative nature of the 
composition law is a consequence of the multiplicative nature of transfer matrices when put 
in series. If the cross section of the conductor is finite the transfer matrix acquires more 
modes. The parallel composition law for classical conductance treats all modes equal and 
so do we. To this end, we first have to define an appropriate conductance for a quasi- ID 
system for which the parallel composition makes sense. The notion of intrinsic resistance R 
and conductance g = Br 1 , respectively, is helpful. In ID we have B = (1 — T)/T which can 
take all positive values. By Eq. ( |165| ) the series composition law reads 

B{bL) = (l + B(L)) b - 1. (166) 

We adopt this law for any series composition in quasi- ID when we put b blocks of size L d 
in series. That this makes sense can be seen by considering the cases with R(L) ^> 1 and 
B(L) <C 1, respectively. For weak resistivity of each block, the resistances behave additive 
like in classical resistor networks. However, for strong resistivity of each block we recover 
the multiplication law leading to strong localization which is consistent with the finding of 
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strong localization in quasi-lD. Once we have put the b blocks in series, forming a quasi-ld 
conductor of cross section L d ~ x and length bL we put b d ~ l of such conductors in parallel, 
resulting in a cubic conductor of total volume (bL) d . For this composition we assume the 
classical parallel composition law to be valid [], 



R(bL) 



■i-d 



(l + R(L) 



1 



(167) 



The application of parallel composition to arbitrary scaling factors can become quite inade- 
quate. However, we will now restrict to very small scaling factors b = l + (,( = 5L/L<^l 
since we are only interested in the change of R under a small increase of system size. By 
this we can easily derive the corresponding /5-function for R, 



d) H — In M + R 



(168) 



dhiL R 

For easier comparison with the scaling theory of Abrahams et al. we reformulate the result 
in terms of the intrinsic conductance 

ding 



(3(\ng) 



-(!+£) In l + j 



(169) 



dlnL 

One sees immediately that this /^-function has the same asymptotic behavior as in the 
scaling theory by Abrahams et al, i.e. for g ^> 1, j3(lng) « d — 2 — g~ 1 /2, and for 3 C 1, 
(3(\ng) »s lng. Furthermore, it has no zero for d < 2 (see Fig. 21) and for d = 2 + e 
with e <C 1 one recovers the results g* = l/e = z/ofa much more involved field theoretic 



calculation ||133|| which is expected to be exact in the limit e < 1. In addition, one finds for 
d — 3 a value of v w 1.65 which is close to known results [y ~ 1.4 |62] , |64]j65[| ) of numerical 
calculations for 3D models . 



C. Scaling Of Distribution Functions 

We are now going to address the flow of distribution functions. As a guiding example we 
consider the ID conductor. Here localization on the scale of the mean free path l e dominates 
and we can expect to get an almost complete description. 

The composition law for transmission, Eq. (0), translates to that of intrinsic resistance, 

R, 

R l2 = 2RxR 2 + 4+^2-2 cos <f>Jk\R~ 2 + (l + R 2 ) (170) 



4 Note, that this procedure is not capable of neither strong magnetic field nor spin-orbit interaction. 
In both of these cases a quantum parallel composition has to be constructed. 
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We can use this to derive an evolution equation for the distribution function V(R, L). Since 
the procedure is of general importance we describe the essential steps (cf. ||112|| ): a stochastic 
quantity X(t) depending on time t may increase in time St by 8X(t) — X(t + St) — X(t). 
If the stochastic process is Markovian, i.e. it is determined by transition probabilities from 
one time step to the next and is not influenced by the history of the process the essential 
information is contained in the moments, ((SX(t)) n ). If, in addition, only the first two 
moments scale linear with St (while all higher moments scale with higher powers) a Fokker- 
Planck equation (FPE) results 

dV(X;t) d[V l {X)V{X-t)\ d 2 [V 2 (X)V(X;t)] 



dt 



OX 



dx< 



(171) 



where V 1 (X) = (5X(t)) /St and V 2 (X) = ((5X(t)) 2 ) /2St (for St -> 0) are called drift and 
diffusion function, respectively. 

In our problem, the stochastic quantity is R(L) and by Eq. ( |170| ) we obtain 



5R(L) = 2R(L)R(SL) + R(SL) - 2 cos <p^R(L)R{SL) (1 + R(L)j (1 + R{SL)) . (172) 

The average is to be taken over random phases and over the strip-resistance R(SL) which 
introduces the mean free path (R(5L)) = SL/l e . This definition is often referred to as local 



weak scattering condition and is the basis of a FPE approach to the statistics. Equation (172) 
shows that only the first and second moments are linear in SL while higher moments scale 
with higher powers. Thus, a FPE can be derived for V(R, L) with drift and diffusion given 
by 

2R + 1 , R 2 + R 



V 1 



The FPE can be brought to the form p5 

dV(R; L) 



dL 



ld_ 

k OR 



[k 2 + r 



dV(R; L) 
OR 



(173) 



(174) 



This distribution depends on only one parameter, the mean free path l e . Two limiting cases 
can be considered now: (i) R 1 and (ii) R ^> 1. 

For case (i) the solution is an exponential distribution 



V(R; L) = yexp (-Rl e /L 



(175) 



with average value (^R(L)J = L/l e . This solution corresponds to system sizes L ^ l e . 
For case (ii) the solution is a log-normal distribution 

.2" 



V{R] L)dR 



1 



AnL/l e 



exp 



InR-L/l, 



2(2L/l e 



dlnR 



(176) 
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with a typical value R t = exp (InRj = exp (L/l e ). This solution corresponds to localization 
with L ^> l e , and a typical localization length of £ = 2/ e , consistent with our previous 
finding (see Eq. fl94|)). Interesting is the fact that the second parameter of the log- normal 
distribution, the log- variance 2L/l e , is simply related to the log- average value 

(\nR- (hii?)) 2 ^ =2(ln^) = 2L/l e . (177) 

The log-normal form of the distribution for L ^> £ is also consistent with the general 
statement, Eq. (|103|) , that the localization length ( ^) corresponds to normally distributed 
quantities that become self averaging in the limit L — > oo. 

The distribution V(R, L) is broad even for L <C £ (since an exponential tail already 
introduces a strong growth of moments oc n!) and causes large fluctuations for R. The 
observation of broad distributions is central to mesoscopic systems and therefore we make 
some general statements about them. Following Shapiro [ |116|| we formulate the concept of 
one-parameter scaling for broad distributions as follows. 

A distribution of a physical quantity X that depends on system size L and a set of 
initial parameters {a n }, V(X; L; {a n }), obeys one-parameter scaling if (for large L) it is 
approximately a function of only X and one scale dependent parameter a^, 

V(X;L;{a n })^F(X;a L ), (178) 

and the flow of the parameter ce L is determined by a /3-function 

That means that the whole set of initial parameters {a n }, needed to specify the initial small- 
scale distribution, disappear into the single parameter ai ({a n }). The exact distribution does 
contain all the information about the initial set of parameters; however, the information can 
be accumulated in the very far tails of the distribution. By very far tail of a distribution 
we denote that range of values X < X\, X > X% that is not statistically significant to some 
small but fixed probability 5 

J dXV(X;L) = 1 -5. (180) 



By Eq. (|178 ) it is meant that the approximation for the main part of the distribution is 



exact up to corrections of order 5. For sufficiently large L it can be taken to be arbitrarily 
small. Still, the moments of the distribution might be dominated by the very far tails of the 
exact distribution. Therefore, one cannot expect that in general the parameter ol is simply 
related to moments of the distribution. 

As a consequence of one-parameter scaling there should exist a universal limiting distri- 
bution at a critical point 
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lim V(X;L) = V*(X) (181) 

which is fixed by the critical value a* = lim ai of the scaling variable. 

A distribution is said to be broad if its moments (X n ) grow like exp/(n) with f(n) 
increasing stronger than linear with n, or if its moments might even diverge. This notion 
refers to the fact that the corresponding distribution of the sum of N independent realizations 
of the variable X will still deviate substantially from a Gaussian, even for very large N. 
According to this definition of broadness a weak case of a broad distribution is realized by 
the exponential distribution (X > 0) 

V{X) = -±rexp{-X/{X)) (182) 

for which the moments grow as (X n ) oc n\ ~ exp(nlnn). The prototype of a broad distri- 
bution is the log-normal distribution 



V(X)dX = exp 

V27TO" 



' (lnX-(lnX)) 2 



2a 2 



dlnX (183) 



the moments of which grow as (X n ) cx exp(o" 2 ra 2 /2). It represents the generic distribution 
for random multiplicative variables. 

After these general remarks we would like to calculate the conductance distribution in 
higher dimensions including the LD transition. However, this is difficult. An interesting 
attack has been undertaken by Shapiro JTTEfl . He used the classical parallel composition law 



to derive Fokker-Planck equations for the conductance distribution in higher dimensions. 
However, the uncontrolled approximations involved overestimate fluctuations. Analytical 
approaches are usually restricted to dimension d = 2 + e where the critical conductance 
is of order e -1 and perturbation theory together with renormalization techniques can be 
used. We do not enter this technically involved subject and we will briefly review the main 
conclusions in Sects. |VIII|JD^ . 



Here, we will turn to local quantities and address the question of power law scaling 
relevant to the delocalized phase and the LD transition. 

In the delocalized phase wavefunctions are extended all over the system. The correlation 
length £ c in this phase can be estimated from 

gfo) = Oil) . (184) 
Taking the metallic limit, Eq. fl29|), one finds that (d > 2) 




(185) 

is a truly microscopic scale. Thus, one can expect a power law scaling of moments of the 
box probability P(E, Lb) for £ c <ti Lb L. In the critical regime of the LD transition the 



51 



correlation length £ c is larger than system size L and one can thus also expect a power law 
scaling of moments of the box probability P(E, L b ) for I <C L b <C L <C £ c (where / stands 
for microscopic scales). Therefore, we will now consider distribution functions the moments 
of which display power law scaling. To simplify notations we consider a random variable 
X G [0, 1] and a scale A G [0, 1] with distribution V[X\ A). We assume that for sufficiently 
small A <C 1 the moments obey power law scaling 

i 

(X q ) x = J dXV[X- \)X q = c q \ f{q) (186) 
o 

reflecting that scales in the system are separated. Here c q and f(q) are treated as functions 
of arbitrary real values q. 

The power law scaling condition is very restrictive and allows for a number of conclusions: 

1. Normalization requires 

c = l, f(0) = 0. (187) 



2. Since f(q) is independent of A one can calculate its derivatives and send A — > 0. By 
this one can conclude two strong inequalities 

r\q) > 0, f"(q) < (188) 



where the prime means differentiation with respect to q. 

3. Introducing the variable Y = — lnX G [0,oo[ one can see that the moments are the 
Laplace transform of the function e~ Y V(X(Y); A), 

oo 

(X*) x = JdY 



Thus, one may reconstruct the distribution from the moments by inverse Laplace 
transformation (if the latter can be carried out). 

4. The case of strictly linear f(q) = ^q (7 = is not interesting) allows for a global 
rescaling of the distribution function by its first moment. 

* P ( X > X ) = Jt i *(§) > Wa = ^, (190) 

and the coefficients c q are determined by the scale independent distribution V(Z) 

1 

c q = J dZV(Z)Z q . (191) 




e- y V(X(Y);X) 



-Yq 



(189) 
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A linear behavior of f seems to be natural for power law scaling; all moments scale 
as powers of the first and only one critical exponent describes the scaling of all of 
them, including that of the distribution function. The rescaled distribution function 
is now scale independent and can take any form, e.g. it can be a broad distribution in 
the sense defined above. Examples are the Rayleigh and Porter-Thomas distributions 
found for the box-probability in the ideal delocalized phase (see Eqs. ( |140| , |141| )). We 
refer to the situation of linear f(q) as the gap scaling situation. 

5. The linear case is not the most general one. For finite curvature in f(q), f"(q) < 0, 
the coefficients c q can be treated as constant since the essential q dependence that 
dominates for small A comes from f(q). We will see in the following that the dis- 
tribution function V(X, A) can be essentially reconstructed from f(q) by a Legendre 
transformation and it is an almost log-normal distribution with special features. One 
can write 

V(X; X)dX = M\-~ s{a) da , a := ^ (192) 

In A 

where the normalization M has only a weak A dependence (e.g. ~ In A). The essential 
information is now contained in the function f(a) which is the Legendre transform of 

%), 

/ (a(q)) = a(q)q - f(q) , a(q) := ^ > . (193) 

Before we outline a derivation of the statement we make further comments. 

6. In contrast to the gap scaling situation with arbitrary form of the scale independent 
distribution we have now the situation that the form of the distribution is fixed by 
a spectrum of scaling exponents. This situation has been phrased (in a bit more 
restricted context) multifractal scaling situation (for a review see | |105| ). The spectrum 
has several analytical properties following from the Legendre transformation. It has a 
unique maximum at (a(0),0) and a slope of 1 at a(l). The derivatives fulfill 

df d 2 f 1 , . 

Close to the typical value 

X t (A) = exp((lnX) A ) oc A a ° , a = a(0) , (195) 

the function f(a) is well approximated by a parabola and V(X; A) is approximated by 
a log-normal distribution 
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V(X; X)dX « A/" exp 



In 2 (X/X t ) 
2 ( x 2 (-lnA) 



dlnX 



(196) 



where a 2 = -l/f"(0) is a positive constant. 



The derivation of Eg. ( |192|) relies on the fact that a finite curvature of r(g) allows for a 



stationary point evaluation of the integral in Eq. ( |186| ) for each moment (X q ) x in the limit 
A — > 0. Writing X = X a and transforming to the variable a, with the distribution V(a, X)da = 
V(X, X)dX, yields the integral relation 



oo 



c q exp [(In A) f(q)] = da exp [— In A (F(a, A) — aq)] (197) 



o 



with 



F(a, A) : 



In (f>(a,X)) 



(198) 



In A 



The function F(a, A) — aq can have a maximum with respect to a (for A < 1) that is only 
weakly dependent on A. It yields a sharp maximum in the integrand (due to the large 
pre-factor -In A). Then, the integrand is determined up to 0(l/(— In A)) by the value of the 
integrand at this maximum and one finds that f(q) is the Legendre transform of F(a,X). 
Consequently, F(a, A) must have a very weak A dependence (e.g. ~ In | In A|). In case that 
F(a, A) — aq has no maximum with respect to a that is only weakly dependent on A the 
corresponding f(q) will not have finite curvature, but the gap scaling situation applies. 

The phrase multifractality has its origin in the interpretation of X(X) as a local measure 
defined for a geometrical object embedded in a d- dimensional cube of linear size L. The scale 
A is given by the ratio (Lb/L) d where is the linear size of a box on which the measure is 
calculated. Due to normalization X(l) = 1 the average of X(X) is the inverse of the number 
of non-empty boxes and thus yields, for A <C 1, the fractal dimension Dp of the geometrical 
object 



^(i/d)(D F +T(q)) w \y eie r (g) defines a set of g-dependent generalized dimension D(q) by 
r(q) =: (q — l)D(q), has led to the notion of multifractality. Especially, one has r(0) = —Dp, 
t(1) = 0. The Legendre transform of r(q) then defines the so-called multifractal /(«)- 
spectrum. It can be interpreted as the fractal dimension of that subset of the object where 
the measure scales locally as (L^/L) 01 (see ||60|| ) . The general analytic properties of multi- 
fractal exponents are visualized in Fig. 22 

We see that, for measures X(X), we have the following identification 




(199) 



The fact that higher moments {X q ) x may scale not simply as oc X <yDF ^ d ^ >q but as oc 
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f(q) = \{D F + r(q)) , f(a) = ~ (-D F + /(a)) , a = - . (200) 

Under quite general conditions, f(a) (or f(a)) is also bounded from below by —Dp/d 
(or t(1) = 0) and terminates with infinite slope at the endpoints of the finite interval 
[a(oo), a(— oo)] (or [a(oo), a(— oo)]). In many cases the parabolic approximation of Eq. (|196|) 



extends to values a > a(l) and the constraint on the slope of f{oi) at a(l) allows to write 
a parabolic approximation (PA) for f(a) that depends only on one-parameter (besides D F ) 

tt'o 

/ (a) = ( a ~ a f r (201) 
4 (a - £>f) 

For example, the box-probability P(E; Lb) considered in the LD-transition problem 
(Eq. (|123|) ) is of this type of variable, since the squared amplitude of the wave function 
defines a local measure. The fractal dimension Dp = d is always trivial. In the idealized 
delocalized phase we found the gap scaling class and r(q) — (q — l)d. The interesting 
information is contained in the shape of the corresponding scale independent distribution. 
For the LD transition, however, it is known from a calculation by Wegner [|133| , p^1 for 2 + e 



dimensions that the multifractal scaling class applies and D(q) ^ d. The PA he calculated 
yields 

a - d = e = l/g* . (202) 

In summary, we found that the restriction to power-law scaling leads to two distinct 
classes of distribution functions. The gap scaling class is characterized by a single scaling 
exponent and a related scale independent distribution of arbitrary shape. The multifractal 
scaling class is characterized by a spectrum of scaling exponents and a shape of the distri- 
bution close to log-normal. The distribution of local box-probability in the ideal delocalized 
phase falls into the gap scaling class. The possibility of the multifractal scaling class is 
interesting for the LD transition point, because it opens a possibility to find an appropriate 
order parameter in terms of the LDOS. If, as it can happen in the multifractal scaling class, 
average and typical DOS scale differently with respect to system size, the inability of the 
average DOS to work as an order parameter must not be shared by the typical DOS. In 
Sect. ILX B| we will see that, indeed, this is the case. 



VIII. STATISTICS IN REALISTIC MODELS 

In this section we will review known facts about the statistics of conductance, energy 
levels and local quantities for the delocalized phase with large typical conductance g 3> 1 
and for the localized phase with finite typical localization length £ 3> l e - The onset of broad 
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distributions with log-normal tails in the delocalized phase is most interesting here. In 2D 
it leads to multifractality of the LDOS. In order to interpret this phenomenon we consider 
the conductance distribution for quasi-lD systems for which the number of channels is large, 
but the width is still of the order of a microscopic scale, e.g. the mean free path. Under 
these conditions a Fokker-Planck approach works that reproduces the UCF phenomenon, 
but the long tails are absent. We come to the conclusion that the log- normal tails are 
a precursor of the LD transition and do only occur if the relevant scaling variable is not 
infinite. To study more realistic quasi-lD systems the method of Finite Size Scaling (FSS) 
will be outlined. In this method a scaling variable is defined which is close to the definition 
of typical conductance. 



A. Onset Of Broad Distributions 

In the delocalized phase of a finite mesoscopic conductor where the average conductance 
is large (g) ^> 1 the Buttiker formula tells that many (O(g)) open modes are occupied. 
Thus, the conductance is a sum of many random numbers % ^ 1. Naive application of 
the central limit theorem suggests that the conductance distribution will be a Gaussian 



centered around (g) with a relative deviation y((5g) 2 )/ (g) ~ l/J(g) where 5g = g — (g). 
However, this reasoning ignores the strong correlations between the transmission eigenvalues 
%. This correlation leads to the occurrence of universal conductance fluctuations (UCF) (cf. 
Sec. [II B|) , i.e. the variance is ((Sg) 2 ) ~ Oil) and does neither depend on the actual value of 
(g) nor on the system size L. Taking UCF into account the conductance distribution seems 
to follow one-parameter scaling: it is Gaussian with scale dependent average value (which 
is also typical here) and a universal variance. 

\2" 



V(g; L) oc exp 



(9 - {9)i 



(203) 



2 ((W> 

This was the common belief until the work of Al'tshuler et al. ||. In this work it was 
shown by a combination of diagrammatic perturbation theory and renormalization group 
techniques that the conductance distribution deviates from a Gaussian for \5g\ \J (g). 
For these tails the distribution crosses over to a log-normal form (see Fig. 23). This was 
concluded from the calculation of high moments which showed the characteristic dependence 

({5g) n ) oc e un2 . (204) 

Here u ~ 0(1/ (g)) indicating that a finite value of conductance is responsible for the effect. 

The calculations were carried out for 2D (with localization lengths much larger than 
system size) and d = 2 + e (in the delocalized phase). In the work by Al'tshuler et al. it 
was conjectured that their finding will penetrate to the LD transition point, i.e. that the 
conductance distribution will become very broad. 
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The observations of these log-normal tails led to questioning the whole idea of one- 
parameter scaling. Later, Shapiro and coworkers [B6J showed that these tails are not in 



conflict with one-parameter scaling theory as defined in Sees. |VII C| . On the basis of the 
2 + e results they concluded that a critical conductance distribution will appear with long 
power-law tails. Such tails lead to divergence of low moments and exclude any one-parameter 
scaling theory in terms of moments of the conductance. 



Later works ||101|, 50,j97 confirmed the calculation of Al'tshuler et al. and showed its 



principle validity for 3D as well. The ID and quasi-lD cases need special attention since 
strong localization dominates. 

In the localized phase the Buttiker formula tells that the conductance is dictated by the 
smallest value z/i corresponding to eigenvalues of M^M (cf. Eq. (|86"D ). This value is known 
to be Gaussian distributed and thus it is no surprise that the conductance distribution, 
for £ <C L, g ~ exp (— £,/L) <C 1, is of the log-normal form. This reflects the Gaussian 
distribution of inverse localization lengths. The typical localization length ((z^) -1 2L) defines 
the typical conductance as the geometric mean g t = exp ((lng)). Here, in principle a second 
parameter besides g t might occur in the distribution, namely the variance of v\. However, 
as happened already in the ID case (cf. Sect. |V11 C| ), all known results indicate that this 
variance is simply related to the average value (cf. e.g. and thus only one-parameter 
seems to describe the distribution accurately, 



V(g; L)dg cx exp 

with 



(ln ^ (ln ^ )2 "rfln,, (205) 



2al 



a\ oc (\ng) L . (206) 
The behavior of the conductance distribution close to the LD transition will be discussed in 

Sect, pen 



Concerning the level statistics we restrict our discussion to the two-level correlation 
function R(s). The energy separation sA can be associated with an inverse time scale. 
The behavior of R{s) is related to the question how a wave packet evolves in time h/(sA) 
T3| , |32| . For times larger than the diffusion time t o the wave packet explores the whole system 



several times and it is not surprising that here the standard random matrix theory results 
are approximately valid. In the standard random matrix theory for the idealized delocalized 
phase R(s) only depends on the level distance s measured in average level spacing units. 
A finite, but large, conductance is related to the diffusion time via the Thouless energy 
-&rh = h/tn 3> A, the relevant transport energy scale. Consequently, the dimensionless 
conductance g (here g can be taken as the average value) sets a scale on the s-line. Indeed, 
one finds different level correlation functions for s <C g and for s ^> g. In the first case, as 
expected, the function R(s) is approximately given by the standard random matrix theory 



57 



result and only small corrections (proportional to powers of 1/g) occur. For example, the 
result for systems with broken time reversal symmetry reads 



. /sin7rs\ Ad . 9/ , /„„„x 

V TTS J 7T z g z 

(Ad a constant depending only on dimension and symmetry). 

The case s 3> g has been investigated by Al'tshuler and Shklovskii M and the finding is 
(for s not exceeding the inverse of the microscopic mean free time r and possible oscillations 
on the scale s ~ 1 are ignored) 

R ^ = g^- d /2 ( 208 ) 

where Cd is a numerical coefficient that depends only on d and on the symmetry class. The 
result can be understood by supposing diffusive motion of wave packets in the system. 

In 2D there exists a broad cross over regime between the two asymptotic regimes of 
Eqs. ( P07| , P08| ) f81fl . This can be understood as follows. The 2D metal resembles a critical 
system in as much as the corresponding /3-function is close to 0. The strong localization that 
will finally dominate comes along with very large localization lengths. It has to be expected 
that the level statistics will change on approaching the LD transition^. A similar behavior 
can be expected to occur in the 2D case. 

The level-correlation function for the localized phase with vanishing localization length 
is given by the Poisson statistics, Eq. ( |119| ), i.e. no-correlation for s ^ 0. To take a finite 
localization length into account which is, however, much smaller than the system size L, 
£ <C L, it is helpful to think of dividing the system into subsystems of volume £ d . The 
level-correlation function R(s;L) will then obey a scaling form [|H| ,|3| 



d 



R{r,L) = \±) F(s(Z/LT). (209) 
In the thermodynamic limit the level correlations vanish, lim R(s, L) = 0, in accordance 

L— >oo 

with the Poisson statistics being the limiting case. 

The scaling form is in agreement with the one-parameter scaling theory since the variable 
£/L is related to the typical conductance in the localized phase, g t = e~ 2L ^. 

The form of the scaling function F(x) can be obtained for x being large or small compared 
to 1. For x ^> 1, the corresponding wave-packet dynamics probes only the region within 
one localization volume and is, thus, diffusive. Consequently one recovers the Al'tshuler- 
Shklovskii behavior 0] 



3 We will discuss some of the consequences in Sect. IXC. 
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F(x) oc x-( 2 -W 



(210) 



For i < 1 but | In x | < L/£ regions much larger than one localization volume are probed 
and the only correlation between levels comes from weak couplings between such volumes 
|| resulting in a logarithmic dependence on x 

F{x) oc (lnx) d . (211) 

We now turn to the statistics of local quantities. Again, for a good metal with g ^> 
1 the distribution of local quantities is essentially given by the standard random matrix 
theory (cf. Eqs. (|140| , |14lD ). However, the long tails are different and the finite value of g 
becomes important. We concentrate on the box-probability P(E; Lb) in the following. This 
generalizes to the reduced LDOS p(E;Lb). We introduce the normalized quantity 

t:= <mt' < p < B v=«/^ < 212 > 

and consider the case g ^> 1. The scale / stands for microscopic box size. The threshold for 
the standard random matrix result turns out to be t ~ l/v^- For t <C l/y/g the Rayleigh 
and Porter-Thomas distributions are valid up to small corrections of 0(l/g) (cf. p8|). 

For t > 1/ ^fg a number of authors [ |101| , [57]| , |9"7| , |120[ | have found the onset of broad distri- 
butions which in 2D and 3D develop long tails of the form 

V{t) ~ exp (-A d \n Bd (t/t d )) . (213) 



Here A d depends on the dimensionality d and on the disorder strength and B 2 — 2. For 
B3 values of 2 and 3 have been predicted. The value td is a typical value serving as ref- 
erence point. The precise expressions are currently under debate (because the authors of 
101| , |50| , |97| , |120|| use different approximative calculation schemes). However, what seems to 



be clear is the occurrence of long-tails of the general form Eq. fl213| ). For the 2D case there is 
also agreement that A2 oc g/(lnL/l) and also td oc g^ 1 ln(L/Z). Due to this scaling behavior 
the result can (for / w l e ) be identified with a multifractal distribution within the parabolic 
approximation (PA, cf. Eq. ( |201| )) 

V(t) - exp f_ ant-(ao-2)ln(//L)) 2 \ 
() P { -ln(L//)4(a -2) ) {2U) 

where the essential parameter ao — d is here (d=2) proportional to the inverse of the con- 
ductance 

a — d oc g^ 1 . (215) 

Thus, states in a good 2D metal (with the localization length being much larger than the 
system size) are mult if r act ally distributed. Since the metallic situation in 2D resembles a 
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critical regime of a LD transition the finding is, after Wegner's pioneering work ||133|| and 



the work by Al'tshuler et al. M a further analytical hint at the multifractality of critical 
eigenstates. 

The case of ID and quasi-lD conductors need special attention. In ID the localization 
length £ is of the order of the mean free path £ ~ Z e and thus there is no room for diffusive be- 
havior. For the quasi-lD case one has to be a bit more precise in fixing the sizes. The length 
L is treated variable and the cross-section Lf -1 is kept fixed. However, the dimensionality 
is in general essential for the determination of physical quantities such as the localization 
length. A strong simplification arises for the case when L t pa l e but a huge number N c 3> 1 
of channels is occupied. Under these conditions the dimensionality becomes inessential and 
the only kinematic quantity besides L is the channel number N c . We refer to this case as 
the idealized quasi-lD case. For the idealized quasi-lD system one can show that 

£(AT C ) = ANJ e (216) 

with A a numerical constant of order 1. We have anticipated this result already in Sect. |V B 



and will sketch a more reliable derivation in the following subsection. The fact that one can 
reach a metallic regime by 

l e < L < £(Ag (217) 

makes the idealized quasi-lD system more interesting than strictly ID. On the other hand, 
the idealized quasi-lD system is still tractable by analytical methods (cf. e.g. ||121|J58|] ). 

In metallic idealized quasi- ID systems one has an average conductance (which is also 
typical) 

g = A'Z/L, (218) 

The distribution of t is, for t < 1/ again given by the standard random matrix theory 
with small corrections of 0(1/ g) (cf. PSfl). However, for t ^> 1/ yfg the situation is different 
as compared to the 2D and 3D case. One finds no log-normal tails but stretched exponential 
tails (at least when the box is shrinking to a point) 

V(t) ~ exp (-B^JcTg^j , B, C ~ 0(1) . (219) 

Note, that g = A£/L is the quasi-lD conductance. 

In the localized regime of idealized quasi-lD systems, ^ C I, where the typical conduc- 
tance is g oc e~ 2L ^, one finds a similar behavior for t 3> (£/L) _1 



P(t)~exp[-BJCt(.£/L)) . (220) 



Note, that £/L is related to the typical conductance here. Both equations, (|219 , |220|) , have 



£/L as the reference point. In order to have t ^> (£/L) 1 one has to have £ ^> l e . Thus, 
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the stretched exponential tail does occur when the box-size is much less the localization 
length. In the case where the box-size becomes much larger than the localization length the 
previously obtained results of the idealized localized phase apply (cf. Sect. |VI A|) . 



The tails of the distribution V(t) that have been observed in idealized quasi-lD (stretched 
exponential), in 2D (log-normal) and in 3D (log-cube) have all been interpreted in terms of 
anomalously localized states (or pre-localized states) P jl01| , pB| . The notion refers to anoma- 
lously rare states with sharp peak in the amplitude on top of an extended backround (in 
3D delocalized phases) or localized states with anomalously short localization radius (in 
d < 2). The question arises what the tails might say in view of the LD transition. The 2D 
case seems to be clear. Here the log-normal tails correspond to multifractal scaling to be 
expected once the correlation length is much larger than system size. This interpretation is 
in accordance with the finding of log-normal tails in d = 2 + e. The distributions become 
multifractal when the conductance is close to its critical value. In 2D the correlation length 
is the localization length which is much larger than the system size in the g ^> 1 regime, and 
thus the states are critical in the sense that the system size is less the correlation length. In 
3D with g ^> 1 the correlation length is (see Eq. ( |185|) ) microscopic and there is no reason to 
expect multifractal scaling. In quasi- ID with g ^> 1 the quasi- ID localization length is much 
larger than the system length L, however due to the strong difference between L t ~ l e and 
L, this length cannot be considered as the correlation length. In strict ID the correlation 
length is the ID localization length which is again microscopic and no critical behavior of 
eigenstates can be expected on scales larger than l e . 

In this context it is worthwhile to discuss the findings about the LDOS for open sys- 
tems. These are defined by surrounding the system with an ideally conducting medium. In 
that case we cannot talk about eigenstates of the system. Still the LDOS as given by the 
imaginary part of the Green's function is a meaningful object. The system itself does not 
conserve total probability which formally breaks hermiticity of the effective Hamiltonian for 
the system. For the 2D and 3D case, however, the results for the distribution of LDOS are 
qualitatively the same as for the distribution of the restricted LDOS in the closed systems. 
A qualitative difference occurs for the quasi-lD case. In the open quasi-lD conductor the 
LDOS shows log-normal tails, 

V(p)~exp(-Agln 2 (p/v)) , (221) 

where v is the average DOS, A is of order 1, and p/v 3> 1- This equation is also valid 
in strict ID. The log-normal tails are not due to the statistics of the eigenstates of the 
corresponding closed system, but can rather be interpreted as reflecting the distribution of 
current-relaxation times for the open system. The latter have a similar log-normal distribu- 
tion. The broadness of this distribution is due to the presence of anomalously long-living 



states |98 



In order to better understand the idealized quasi-lD system we will focus in the following 
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subsection on the description of transport in idealized quasi- ID systems. One particular 
question we like to address [25|: is it possible to obtain similar log-normal tails for the 
conductance distribution in the metallic regime as those found in the d = 2 + e case by 
APtshuler et al. §? 



B. Statistics in Quasi-One-Dimension 

The transport properties of a quasi- ID system can be obtained within a transfer matrix 



modeling as outlined in Sect. [V A| . The characteristic parameters are the length L, the 
number of channels N c related to the cross-section Lf~ x and the mean free path l e . 

In close analogy to the procedure that led to a FPE for the intrinsic resistance in ID 
(Eq. ( |174| ) one can derive a FPE for the whole transfer matrix M. The essential assumptions 
for this to be possible are: 

1. The locally weak scattering condition 

l e sl^o 5L K ' 

2. The statistical independence of strip transfer matrices defined for adjacent strips. 

Fixing the statistical properties of a single strip transfer matrix yields a FPE for V(M; L). 
One particular model turned out to be analytically tractable to a large extent, the so- 
called isotropic quasi-lD model. In this model |^2| it is assumed that forward and backward 



scattering within a strip S-matrix is equal for all channels. This crucial simplification means 
that there is no diffusion in transversal direction. The modeling thus restricts the width L t 



to be less than the mean free path l e ||93|| . Therefore, an idealized quasi-ld model in the sense 
of the foregoing section results. Furthermore, the most explicit results can be obtained in 
the limit of large channel numbers, or more precisely when N c — > oo but N c l e /L remaining 
fixed and finite. 

In terms of the polar parameterization Eq . (|87D the isotropy assumption is equivalent to 
the assumption that the unitary matrices u % are distributed isotropically within the uni- 
tary group. Similar to the standard random matrix theory for Hamiltonians the statistics 
of eigenvectors is treated as being trivial and decoupled from the statistics of the radial 
parameters \. 

A consequence of the isotropy assumption is the fact that a closed FPE can be found for 
the distribution of the radial parameters \ alone. These radial parameters determine the 
conductance (Eq. (|86|)) 

* = E* = £r^-. (223) 
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The derivation follows the same principle steps as outlined in Sect. |VII C| and the final 
result is the well-known DMPK-equation E2 



dL 7 Y d\i d\i 

J= II |A m -A n |^, 7 = /3iV c + 2-/3. (224) 

ra<m 

Here (3 = 1,2,4 characterizes different symmetry classes (systems with time reversal sym- 
metry (3 = 1, systems with broken time reversal symmetry (3 = 2, systems with time reversal 
symmetry and spin-orbit interaction (3 = 4) and J is the Jacobian for the transformation 
from the Cartesian parameterization M^i to the polar parameterization of the transfer ma- 
trix. 

The DMPK equation has been solved exactly by Beenakker and Rejaei [[17]] for (3 = 2 



and for (3 = 1,4 formal solutions exist f2"3fl . Common to the solutions is that they represent 



a Gibbs-ensemble of the form already encountered in the standard random matrix theory 
for energy levels, Eqs. (|133| , |134|) , 



VHK}) = exp(-/9W({A 4 })) , 

n ({A,}) = l -Y. U (*n, A m ) + £ V(\) ■ (225) 

m^n i 

For (3 = 2 the explicit form of the two- and one-body potentials read 

1 1 
U(x, y) = —— ln(|x — y\) — — InQarsinh 2 ,/?/ — arsinh 2 -^/^]) , 

V(x) = ^arsinh 2 ^ + 0{l/N c ) . (226) 

It is worth mentioning that, for \ (open channels % « 1), the two-body potential 
reduces to 

U{x,y) = -\n{\x-y\) (227) 

which has been used to model quasi-ld wires, before the exact solution was found (see 
e.g. ||121|| ). A general feature of the Gibbs-ensembles is that physical parameters, such as 
N c , l e and L do only occur in the one-body potential V(x) in the combination N c l e /2L. In 
contrast, the two-body potential is a universal function reflecting universal level repulsion. 

Without solving the DMPK equation explicitly one can study the resulting equations for 
the moments (g n ) by integrating on both sides of the DMPK equation over + Ai)^ 1 )™. 

A perturbative analysis as a series expansion in powers of the single parameter N c l e /L is 



possible in the limit N c l e 3> L ^> / e , i.e. in the diffusive limit (Mello et al. in Ref. [p2|j). 
This yields the result for the average conductance 



(g) = NJJL + C p + ((NJJL)- 1 ) . (228) 
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The constant Cp contains corrections to the classical result (for (3=1 these are negative 
weak localization corrections). Thus, the leading term is the one we have already anticipated 
in Eq. fl95|) . The leading term for the variance ((5g) 2 ) turns out to be universal 



W) = ^ (229) 

and proofs the UCF effect for isotropic quasi-lD. The effort to derive results for higher 
moments by the perturbative method increases drastically. 

In the limit L ^> N c l e it is advantageous to transform to the parameters v% (Eq. (|86f) ) 
since they become the Lyapunov exponents in the thermodynamic limit. Their interaction 
can then be neglected and one ends up with an effective joint probability density of the form 



N 



(/3N C + 2 - p% ( L{(3N C + l-(3i)V 



2L \ l((3N c + 2-(3 



(230) 



demonstrating again the self-averaging feature of Lyapunov exponents. In that limit the 
conductance is g = 2e~ Ul since all higher z/, contribute much less. The conductance distri- 
bution is then found to be of the log-normal form. The typical value and the log-variance 
are related by (N c > 1) Q 

((5\ng) 2 ) = -2hxg t = -2 (lug) = ^ . (231) 
The localization length is then found to be 

f = f3N c l e . (232) 

Before we look at the conductance distribution in more detail we report on the findings 
concerning the distribution of transmittances T a p = |t Q/ 3| 2 and T a = 

They are related to the eigenvalues % = (1 + A;) -1 of tH by unitary matrices u and v 



T otP = u akU* al \J%Tiv k pv* w , T a = J2 \uia\ 2 % ■ (233) 



kl 



These equations point out that transmittances are analogs of wave function amplitudes 
I V^ct (^) | 2 — |£4a| 2 - However, the transmittances do depend also explicitly on the eigenvalues 
%. The statistics of the unitary matrices was assumed to be isotropic in the DMPK approach 
and the eigenvalues follow the DMPK equation. 



6 Notice, that the so-called global transfer matrix approach | 108|J121 yields a factor 1 between 
variance and average value instead of a factor 2 occuring in Eq. ( p31[ ). 
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It turned out that the distribution of transmittances can be calculated exactly |86 . 



Actually, only V(T n ) needs to be calculated in the large N c limit. In a work by Kogan and 
Kaveh Jr|| it was shown that, for isotropic scattering in the large N c limit, the distributions 
of T a and T Q/ g are connected by 

oo 

V(T af3 ) = J dT a T- 1 exp {-T aP /T a ) V(T a ) . (234) 
o 

In the metallic regime NJ e ^> L, (g) ^> 1, the distribution of T a p is of the Rayleigh type for 
the bulk of the distribution (j3 = 2) and develops stretched exponential tails (first obtained 
by Nieuwenhuizen and van Rossum ||102|| ) 

V(T aP ) ~ exp (-2^/55) , T aP > (N c l e /L) 2 . (235) 

It is therefore completely analogous to the distribution of wave function amplitudes in closed 
quasi- ID systems (see Eq. fl219|) ). The corresponding distribution of T a is Gaussian in the 
bulk and develops exponential tails 

V{T a ) ~ exp (-T a ) , T a > N c l e /L . (236) 

These results depend heavily on the isotropy assumption and they indicate that it is very 
unlikely that the distribution of conductance, g = J2 a T a , in the metallic regime can develop 
log- normal tails in isotropic quasi- ID. However, so far the extension to conductance is not 
easy since the correlation between eigenvalues % have to be taken into account. 

In the localized regime L ^> N c l e — the statistics is dominated again by the smallest 
value V\ of i.e. by the largest value T\ of %. This value is log-normally distributed 
(see Eq. ( |231|) ) and accordingly the transmittances are log-normally distributed, too. 

Since the conductance is a linear statistics of the transfer matrix radial coordinates \ 

9 = J2f( X i)> f(x) = (l + x)- 1 (237) 

i 

the problem of a complete statistical description of g reduces to the problem of a linear 
statistics X = J2if( x i) i n a Gibbs-ensemble of a classical gas determined by a universal 
two-body potential (which can be chosen to be symmetric) and a parameter-dependent 
one-body (confining) potential 

V(X)= J d N ^xV({x t })5^X -^fix^ , V({x t }) = Z' 1 exp-(3H({x l }) 
n {{xi}) =^E U ^ x m ) + E Vfa) . (238) 

m,n i 



Beenakker pointed out [IB that the one-body potential can be viewed as a source term 
in the partition sum 
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Z[V}= J d Nc xexp[-[3H({xi})} . (239) 
All cumulants [] of the level-density, p(x) := J2i &{ x 

be obtained by functional 

derivatives 

(^■■■« = iF^F^M- (240) 

Cumulants of a linear statistics X = J2i f( x i) are given by integration 

((X k )) = Jdxt... dx k ((p( Xl ) . . . p(x k )}} /(an) . . . f(x k ) . (241) 
The whole distribution V(X) can be obtained from a modified partition sum Z(k), 

V{X) = J dKe iKX Z(K), (242) 
where Z{k) follows from Z through a simple shift in the one-body potential 

V K (x) :=V(x)+£f(x). (243) 

In order to work with functional derivatives one has to know Z (or Z(k)) as an ex- 
plicit functional of the one-body-potential. Alternatively, the knowledge of the average level 
density u(x) := (p(x)) as a functional of the one-body-potential would be enough. 

Such representation of the average level-density v(x) was developed by Dyson p|. Here 



we outline a derivation |29[ which allows for a systematic check of the range of validity of 
the results. 

We start from the partition sum Z, or Z(k) and proceed by the following steps: 1. The 
Hamiltonian is expressed by the level-density p(x), replacing sums by integration. Here the 
absence of self-interaction is ignored, i.e U(x n , x m ) is replaced by / dx dy p(x)p(y)U(x, y). 

2. With the help of a 5-functional we introduce a field <p(x) that takes the role of the level- 
density. 3. We represent the ^-functional by its Fourier representation which introduces 
a conjugate field ip(x). 4. Now the original integration over the set {x{\ can be carried 
out leaving a field theoretical partition function in terms of two field degrees of freedom, 
4>(x) and ip(x). 5. Due to the two-body character of the original 7i the field 4>(x) can be 
integrated out by a Gaussian integration. 



7 Cumulants ((X n )} are linear combinations of moments of order k < n. While moments can 
be generated from a partition sum Z, the corresponding cumulants are generated by InZ. The 
Gaussian distribution is characterized by vanishing cumulants for n > 3 
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The final result is a path integral representation of the partition sum Z or of the complete 
distribution function V(X), where the integration runs over field configurations of tp(x). We 
concentrate on V(X) which reads 



V(X) = jDMexp-SMX] 



S *] =™ ty| K\ip) + ^ (f\ K \f) Q\i,- X) + F Nc [V + pV\ 



Here we use a scalar product notation, 

(f\A\g) := f dxdyf(x)A(x,y)g(y), {f\g) := 



dxf(x)g(x) 



(244) 



(245) 



K denotes the inverse operator of U (U (x, y) = (x\ U \y)), the functional Q{ijj; X) is defined 

as 



Q(^X) :-. 



P- 1 (f\K\f) 



(246) 



and the functional F^ c is a free energy of N c independent particles with one-body-potential 



F Nc [V> + 0V] := -N c In J dx exp [- (^(x) + /3V(x))) 



(247) 



The omission of the self-interaction can be cured on an effective potential level. By 
shifting the one-body-potential V(x) to V(x) = V(x) — ~U(x,x + A(x)) where A(x) is the 
average level spacing at x, the theory is able to account for those quantities that are smooth 
on the scale of A(x). For example, with the logarithmic interaction, U(x, y) = — In | x — y \ 
the effective potential reads V(x) = V(x) — ~lnu(x). In the following, we will denote V by 
V. 

The path integral in Eq. (|244j) cannot be calculated exactly. However, for large A^ c one 
can use the method of stationary point. To be more explicit we will now use the conductance 
g as linear statistics variable and recall that V(A) contains the pre-factor g = NJ e /L which 
we assume to be large: go ^> 1 since we are looking for the conductance distribution in the 
metallic regime. 



Introducing the mean-field level density as 



N 



i>° + 13V 



9 ' 



d\ exp [- (>°(A) + V(X) 



(248) 



corresponding to the stationary point, 



5S_ 

5^ 



0, 



(249) 
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the mean-field equation rea 



v. 



-K 



V + f3- l \nv°\+f3- l Q{g)K\f) + 



N C K 


1) 


(1 


K 


1) 



Here the kernel K is defined as 

K — K — 

and Q(g) as 



K\ 


1)(1 


\K 


(1 


\K\ 


1) 



Q(g) 



K\l) 



9~9 g 



[l\K = 



p-Hf\K\f) 



with 



g q :=-(f\K F + /TW) + 



N c (f\K\l) 
(1\K\1) 



(250) 



(251) 



(252) 



(253) 



In deriving these equations it has been used that the mean-field level density is normalized 



(l | Vg^J = N c and yields the current g as an expectation value 
Now one can draw the following conclusions: 



9- 



1. For |V(A)| ^> (3 1 lni/° the expression g g equals the average value of g, independently 
of current g. Since V contains the large factor g the inequality is satisfied as long as 



N < iff) , (9) > 1. 



(254) 



2. The stationary point of S is then given by 



S 



1 (9-(g)f 



2/3- 1 (f\K\f) 
where is independent of current g. 



+ S, 



(9) 



(255) 



3. One can also analyze fluctuations around the stationary point and show that they give 
sub-leading contributions to the path integral. Finally, one arrives at the conclusion, 
that the distribution of conductance is, for \8g\ <^ (g), 



V(g; L) = const, exp 
up to 0(lng) corrections in the exponent 



(9 ~ <<7>j 



L 2/3- 1 (f\K\f) 



(256) 



8 This mean- field equation is in agreement with the one obtained by Dyson [43 and Beenakker 
[16] for the special two-body-potentials used in that works. 
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Therefore, our conclusion is that the conductance distribution in isotropic quasi-lD does 
not show long-tails in the regime \5g\ <C (g), but is Gaussian with universal variance. The 
universality is due to the universal two-body potential U(x, y) which reflects essentially the 
universal level repulsion expressed by the Jacobian in Eq. ( |224|) . Recall, that the occurrence 
of log- normal tails in d = 2 + e happened already for \5g\ ^> J (g). 

The absence of log-normal tails in the conductance distribution for isotropic quasi- ID 
can be interpreted as being caused by the isotropy assumption. Relaxing from the isotropy 
assumption the resulting FPE for the radial parameters Aj is not closed. Rather the diffusion 
function depends on the statistics of eigenvectors of M^M. The latter is presently unknown, 



but one can conclude in general |30[ that long tails in the distribution of the eigenvectors 
will lead to long-tails in the distribution of g. In addition, one can conclude that in such 
case the probability density of the radial parameters is no longer given by a Gibbs-ensemble 
with only one-body and two-body potentials. Also the results for transmittances indicate 
that the isotropy assumption expels any log-normal tails for the conductance distribution. 

On the other hand, the isotropy assumption can only be justified for transversal width 
L t being less the mean free path l e . It is instructive to look at the correlation length £ c 
relevant to such idealized quasi- ID systems. Recall, that the correlation length, in the spirit 
of one-parameter scaling theory, is defined as that fictitious system size £ c for which the 
conductance of a rf-dimensional cubic system is of order 1. To attach a correlation length 
to the quasi- ID system requires a meaningful definition of a conductance in a cubic system. 
The easiest way is to use a parallel composition law and define 

g C ubc(L) = g qu£b (L,L t ) {L/Ltf^ 1 (257) 

where g qua (L) is the conductance of the quasi-lD system with length L and cross-section 
L t ~ l . For the isotropic quasi-lD system with g qua (L) ta N c l e /L and L t = Al e , A < 1, we 
thus have 



L 



-2 



g cuhc (L) ~ N c (^-J (258) 

which is extremely large. Consequently, the correlation length is 

£ c « i e iV- 1 /(^-2) . (259) 

Thus, even for d = 2 the isotropic quasi-lD model corresponds to microscopic correlation 
lengths. One should note that it is not the quasi-lD localization length £ = f3NJ e which 
plays the role of the correlation length, but £ c . Both coincide only for strict ID where 
N c = l. 

We interpret the occurrence of log-normal tails in the conductance distribution as pre- 
cursors of the LD transition. For d = 2 + e they are controlled by the one-parameter scaling 
variable g t . With this interpretation the absence of such tails in isotropic quasi-lD is no 
surprise since the appropriate scaling variable oc N c (L/l e ) d ~ 2 is sent to infinity. 
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C. Finite Size Scaling 



Giving up the isotropy assumption allows to study more realistic systems by a quasi- 
1D set-up. True dimensionality effects should enter the physical properties. An analytical 
treatment becomes much more difficult, but a numerical treatment turns out to be very 
efficient. The Green's function can be studied by recursive methods |Kj, putting strips of 



cross-section Lf~ l together. Alternatively, strip transfer matrices can be multiplied ||107|| . 
That quantity which can be obtained with high precision (due to its self-averaging property) 
is the quasi-lD localization length £,(L t ). By studying £,(L t ) as a function of L t one should 
be able to study the localization behavior in a large (cubic) d- dimensional system. 

Our first task is to construct a suitable scaling variable which can be derived from £,{L t ). 
When we increase Lt and the system flows towards the insulating (localized) phase, £(Z^) 
will stop growing as soon as L t is comparable with the (finite!) localization length of the 
infinite d- dimensional system: 

£(oo) = lim i{Lt) ■ (260) 

On the other hand, £,{L t ) will grow forever, if the system floats to the metallic (delocalized) 
phase. 

At the critical point, however, we expect the system to be scale invariant (diverging 
correlation length), i.e. £,{L t ) should scale exactly as L t as soon as finite size effects have 
died out. So one tries 

A(L t ) := i{L t )/L t (261) 

as a scaling variable. 

A further reason for this choice comes from the interpretation that, in the metallic regime, 
an appropriate conductance for the corresponding cubic system is given by Eq. ( [^57| ) where 
# qua (L; L t ) « £{L t )/L yielding 

g cuhc (L = L t ) « A{L t ) (262) 

In the localized regime, where £(Lt) approaches £(oo), A(L t ) represents the inverse logarithm 
of the corresponding conductance 

-\ng cuhc (L = L t )^A- 1 (L t ) (263) 

Now, one-parameter scaling means that away from the critical point (where A(L t ) = A* is a 
constant) A(L t ) does not depend on the correlation length £ c (which can be identified with 
£(oo) in the localized phase) and L t separately, but only on their dimensionless ratio 

(264) 
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Here / is a called the scaling function. In the presence of LD transitions, / is a two-valued 
function. Close to a critical point (where linearizing of the corresponding /5-function is 
appropriate) the scaling function is of the form (cf. Eq. (|160 )) 



f(x) = A* ± A (±) x" (265) 

where the positive (negative) sign marks the delocalized (localized) branch of /: /3(ln A) > 
03(lnA) < 0). 

The general procedure to extract the scaling function from a large number of calculated 
localization lengths £(L t ) is to find a quantity £ c independent of L t such that the entire data 
set £(Lt] t) (where r is a system parameter related to e.g. the disorder or Fermi energy with 
critical value r* = 0) collapses onto a single curve. This curve then defines /. In practice, 
this can be done by plotting the lnA(L t ;r) curves on transparencies and moving them 
around by hand until they fall on top of each other or by a least squares fitting procedure. 
A typical example for a scaling function is shown in Fig. 24. After determining A* one can 
determine the critical exponent v from a fit to Eq. (J265) . 



It may, however, turn out to be difficult to find the critical value A*. For example, it 
can happen that the curves A(L t ; r), for different choices of L t , do not intersect at one point 
r = 0. This can be due to the fact that the system size is not large enough to reach the 
asymptotic scaling regime, i.e. it can be a finite size effect. Assume that, at the true critical 
point t = t* = 0, the scaling variable still shows a dependence on L t of the form 

A(L t ; r = 0) - A* = BL^ Vm (266) 

Here y- m > is a critical exponent. More convenient is BL~[ Vi " = B (L t /^ TT ) ~ V ' IT with B of 
order 1 and an explicit scale £i rr called irrelevant length. Since stays finite the correction 
to A* decreases when L t becomes larger than £ irr . 

lim A(L t ; r = 0) = A* . (267) 

Lt^oo 

In this sense, the length £i rr becomes irrelevant in the large L t limit. However, in a finite 
size calculation one can never be sure that microscopic scales are already much smaller than 
L t and one has to face the fact that corrections to the ideal scaling, Eq. ( |265| ), can occur. 
They can then be analyzed by means of irrelevant length scales, Eq. (|267|) . 

The best studied system is the quantum Hall system in 2D. Although no metallic phase 
exists, the position of the critical point on the energy scale is precisely known for certain 
models due to spectral symmetries. This fact and the two- dimensionality enables to get 
results for critical quantities with high precision (see e.g. f6~8|j72|| ). Also a 3D analog of 
the quantum Hall system has been studied recently pi |. Further results are known for 2D 



systems with spin-orbit scattering (see e.g. ||5]|)T],[IT^]) and for a number of 3D systems (see 
e.g. BUS). 
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IX. STATISTICS AT THE TRANSITION 



The finite size scaling method made it possible to analyze a self- aver aging quantity 
serving as scaling variable in the sense of one-parameter scaling theory. It demonstrates 
not only that one-parameter scaling does work, but also allows for explicit calculations of 
the critical exponent of the correlation length. However, not very much information comes 
out of this for the statistical properties of non-self-averaging quantities like conductance and 
LDOS. At least, by Eqs. (|262|J263|) one can expect that the scaling variable A of the finite 
size scaling method is related to the typical conductance gt- The latter can be identified as 
the mean value in the delocalized phase and as the geometric mean in the localized phase. 
Since, in the delocalized phase the geometric mean is very close to the mean (the distribution 
is Gaussian with yj ((5g) 2 ) / (g) <C 1) it is reasonable to expect that the geometric 

mean 

<7 t :=exp((ln<7)) (268) 

can serve as a typical conductance which in turn is a scaling variable. 

In this section we will firstly consider the statistics of critical eigenstates. After the 
findings of the previous section we already expect them to show the multifractal scaling 
property. We will discuss the implications of this for the LDOS. Several arguments will 
point out that the typical conductance will determine the distributions of local quantities. 
After a brief discussion of critical energy level statistics we collect results for the critical 
conductance distribution. 



A. Multifractality Of Critical Eigenstates 

In finite size systems the correlation length £ c of the electronic states is larger than the 
system size L for a certain parameter range, At, around the critical parameter value r = 0. 
These states are called critical states. In the thermodynamic limit At oc L~ x I v where v is 
the critical exponent of £ c . 



After the pioneering works by Wegner ||133|| and Aoki |0| it became clear that the 



critical wave functions have a multifractal structure (for a review see |7T[] and references 
therein). The entire distribution of local amplitudes and its scaling behavior is encoded in 
the multifractal f{a) spectrum, as outlined in Sect. |VII Cj . The distribution V(P; L b /L) 



is broad on all length scales and close to a log-normal distribution. The most important 
quantity is the maximum position, ao, of f(a). It describes the scaling behavior of the 
geometric mean of what serves as a typical amplitude of a critical wave function. 

In early works Aoki |TTJ gave a nice argument for the multifractal behavior of critical 
wave functions (although at that time the phrase multifractality was not yet common). His 
argument goes as follows. Consider the inverse participation number defined by 
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V 



Jd d r\iP(v)\ 4 (269) 



where Q denotes a d-dimensional region with linear size L. If the wave function if>(r) is 
uniformly distributed - like in the delocalized phase - then V oc L~ d and the participation 
ratio p = {VL d )^ 1 is constant. In the localized regime V ~ and p vanishes in the 
thermodynamic limit. At the transition point where the wave function is extended the 
participation ratio still has to vanish in the thermodynamic limit if the LD phenomenon is 
similar to a second order phase transition for which the participation ratio acts as an order 
parameter. Consequently, V scales with a power d* < d. Wegner had already calculated 
the whole spectrum of exponents for generalized inverse participation numbers within the 
non-linear sigma- model ||133|| . This spectrum was interpreted as a multifractal spectrum 



by Castellani and Peliti [24]]. After extensive numerical work (e.g. [[L09| , |117|j66|J118|1 ) the 



following description of the statistics of critical wave functions ^>(r) is now established: 

At the LD transition the distribution function V{P\L b j ' L) gives rise to the power law 
scaling for the moments, 

([P(L b )Y) L oc (L b /L) d+ ^) , (270) 

where d + r(q) is a non-linear function of q. This non-linearity is a direct consequence of 
Aoki's observation that d + r(2) = d + d* ^ d + d. In practice it turns out that, to a 
good accuracy, the disorder average can often be substituted by the spatial average over one 
wave function for a given configuration. Within numerical accuracy the resulting spectra 
are identical. 

The corresponding (universal) distribution function can be described in terms of the 
single-humped, positive f(a) spectrum, 

V{P; L b /L) dP oc {L b /L) d ~ fia) da , (271) 



where a := InP/ ln(L b /L); As shown in Sect. VII C| f(a) is related to r(q) by a Legendre 
transformation 

f(a(q)) = a(q)q-T(q) , a(q) = dr(q)/dq . (272) 

The parabolic approximation (PA), Eq. (|201|) , contains «o as the only parameter be- 
sides d. This is due to the assumed validity of the PA at least up to \q\ < 1. Equa- 
tion (|201|) corresponds to a log-normal distribution centered around the typical value 
P t = exp < InP >oc (L b /L) a ° with log- variance proportional to — d. A simple one- 
parameter approximation for f(a) which takes into account that the support [a(oo), ot(— oo)] 
of f(u) is finite, is the semi-elliptic approximation (SEA) 



f(a) ss d, 



(273) 

\ etc — d 2 
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To demonstrate that the distribution of local amplitudes of critical eigenstates is encoded 
in f(a) we discuss numerical results for a quantum Hall system (QHS) |]11U[| . The quantum 
Hall system is the most convenient to study critical eigenstates since one knows the critical 
point exactly and it is only two-dimensional which keeps the numerical effort low. 

The wave functions are calculated for the model of independent (spin-less) electrons 
subject to strong magnetic field and disorder. The disorder was implemented by a set of 
5-impurities with random positions and random strengths symmetric around zero. The 
microscopic length scale of the problem is the magnetic length Ib defined by the size of 
an area penetrated by a single flux-quantum, 2ttI b . The representing Hamiltonian matrix 
was worked out in the Landau representation diagonalizing the clean system (no disorder). 
The clean case has the famous Landau level spectrum of highly degenerate and equally 
spaced levels at discrete energy values. The spacing is given by the cyclotron frequency 
Tiuj c = heB/m. Disorder (symmetric around zero energy) broadens the Landau levels to 
Landau bands of non-degenerate states. For strong magnetic fields the broadened Landau 
bands are still separeted, the band- width T being smaller than Tiuj c (see Fig. 25). 

Restricting to one Landau band (it is most convenient to take the lowest) the LD tran- 
sition takes place precisely at the center of the symmetric Landau band (cf. |f75| , fl68|). A 
system of area L 2 is represented by a finite matrix of dimension N = L 2 / (2-kI b ). The nu- 
merical diagonalization yields the eigenvalues and eigenstates for any desired energy window 
within the lowest Landau band. 

In Fig. 26 the squared amplitudes of a wave function from the center of the Landau 
band are shown together with the f(ot) spectrum calculated from these amplitudes. The 
corresponding histogram of the logarithm of amplitudes (measured on a box of size Al 2 B ) 
is displayed in Fig. 27 together with the distribution function calculated from the f{ce) 
spectrum using Eq. ( |271| ). These figures demonstrate that the distribution of amplitudes is 
(i) encoded in the f(a) spectrum and (ii) is close to a log-normal distribution characterized 
by one critical exponent cto = 2.28 ± 0.03 (the average over 130 critical states). 

Figures 26, 27 correspond to a model of 5- impurities resulting in a short range poten- 
tial correlation. Similar calculations for a long range potential correlation (the potential 
correlation length is much larger than the magnetic length) have been performed within 
the Chalker-Coddington network model [f74| |. A corresponding critical wave function and its 



histogram are shown in Figs. 28, 29. They demonstrate the universality of the multifractal 
properties with respect to the potential correlation length. 

To study the spatial correlations of amplitudes for a fixed energy consider the 
q— dependent correlations 

M^(r,L b ,L) := {[P^L^P^L^) L (274) 

where the average is to be taken over all pairs of boxes with fixed distance r = sL b . 

For critical states where the microscopic scale (Ib in our case) and the macroscopic scale 
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(the localization length £ in our case) are separated, one can expect power law behavior of 
in the regime 



l B < L b , r, L < £ . (275) 

Usually, in critical phenomena one studies correlations for infinite system size (and L b being 
microscopic) as a function of r alone. This is justified if, for large enough system sizes 
L, the correlation function is independent of L (for simplicity we neglect any trivial in- 
dependence due to pre-factors in the definition of the observable P). However, this is not 
true in the multifractal case. Multifractality reflects broadness of the distribution function 
V(P, Lb/ ' L) on all length scales. The local box observable Pi(L b ) depends on a large number 
of conditions for the entire system of linear size L, simultaneously. In the context of the LD 
transition coherence at zero temperature is due to quantum mechanical phase coherence of 
the electron's wave function, and disorder introduces a huge number of parameters, e.g. the 
position of point-scatterers. In the multifractal scaling case one has to face the fact that 
AfW depends non-trivially on L, even for L — > oo. We have incorporated such behavior 
already in our general discussion on correlation functions in Sect. [VII A . 



Therefore, we consider the regime Is "C L b < r < L <C £ c and make the ansatz 

M [q] (r, L b , L) oc L x b 2[q) L- y2ig) r- z{q) . (276) 

The task is now to find the scaling relations between the set of exponents x 2 (q), 1/2(1), z(q) 
and the r(q) function. These scaling relations can be derived by requiring consistency with 
the limiting situations (i) where r is of the order of L b and (ii) where r is of the order of L. 



We find [25,71 



y 2 {q) = d + r{2q) (277) 
x 2 (q) = 2d + 2r(q) (278) 
z(q) = d + 2r(q) - r(2q) . (279) 

It is worth mentioning that the sum x 2 {q) — y 2 {q) — z(q) vanishes due to the normalization 
of the wave function. 

The analytic behavior of z(q) according to Eq. ( |279j ) is shown in Fig. 30. In general, it 
is non-negative and asymptotically bounded by the dimension d . To check on the validity 
of Eq. (|279D Pracz et al. took 100 critical states of a system with L = 200Z# and calculated 
M^(L b ,r,L) with fixed values L b = Z B ,4Z B ; L = 200l B . Within the errors the validity of 
the scaling relations in Eqs. (|279| , 278|) could be confirmed. 



To include aspects of local energy statistics as well, let us consider the g-dependent 
correlation of box probabilities corresponding to two different eigenstates with energies E 
and E + u> 

M®(r,L b ,L) := ([Pi(E; L b )m +S (E + w; L b )]«) L . (280) 
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To understand the correlation behavior of non-localized states with respect to the energy 
separation one has to compare the relevant energy scales of the problem. These are the 
average level spacing A and the (frequency dependent) Thouless energy energy E Th (u) 
corresponding to the time a wave packet (formed from states within an energy window of 
width to) needs to diffuse through the system, L 2 = (h/ 'Ej^(uj))D(uj). Here D{uj) is the 
corresponding diffusion constant. According to Chalker |28| , these scales give rise to the 



definition of two length scales depending on the energy separation uo: 

Luj := (uj/E^y^L (281) 
L u :={u/A)-V d L. (282) 

The first length scale, L^, is the typical distance a wave packet will travel diffusively in 
a time h/uj. From this it is natural to assume that correlations between Pi(E; L^) and 
Pi +S (E + lu] Lb) will be present at least for distances r <^ whereas for larger distances 
the amplitudes are uncorrelated. Such uncorrelated behavior is typical for the standard 
random matrix theory approach to extended states in random systems. Recall, that in 
this theory it is assumed that the unitary matrices that diagonalize the Hamiltonian are 
distributed uniformly in the unitary group and no correlations (apart from the unitarity 
property) between different matrix elements occur. Thus, the presence of correlations here 
is an explicit breakdown of the no-preferential basis assumption. In electron systems with 
spatial disorder a preference to some basis is always given. This preference cannot be seen 
in M$ for distances r ^> L u . 

The second length scale, L u , is the linear size of a system with level spacing lu. Two wave 
functions with energetic separation smaller than the level spacing show a spatial correlation 
behavior of its amplitudes similar to that corresponding to one of those wave functions, i.e. 
they are statistically indistinguishable. 

At the critical point of the LD transition the typical conductance becomes a size inde- 
pendent quantity, g*, and with the help of the Einstein relation between conductivity and 
diffusion one finds L u = (g*) 1 ^ d L UJ p8| . Since g* is of 0(1) the two length scales coincide at 
the LD transition. Therefore, one can focus on the role of L^. 

Asking for the scaling properties of M$ in the regime < r < L w < L we make the 
ansatz (cf. Eq. ( |276| )) 

M$ oc L ^ 2(q) r- Z ^L Z ^L- Y ^ . (283) 

Here we have already anticipated that the exponent with respect to r is z(q), as given before. 
As in the case of for a fixed energy by considering limiting situations one finds 

X 2 (q) = 2d + 2T(q) = Y 2 (q) (284) 
z(q) = Z(q) =d + 2r(q) - r(2q) . (285) 
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Now the following conclusions can be drawn: 1. The result for z(q) is the same as in the 
case of zero energy separation. 2. The energy separation u is not an independent scaling 
parameter but appears only in the combination L u /r. 3. The exponent corresponding to 
the box size, 22 = ^(q), is not affected by a finite energy separation but the exponent 
corresponding to the system size L (which is 1/2(9) for zero energy separation) splits up into 
the exponents Z(q) (corresponding to L u ) and Y 2 (q) (corresponding to L for finite energy 
separation). 

In the work by Pracz et al. it was demonstrated that the scaling relations are consistent 
with numerical results. Similar findings have been obtained by Metzler [96] within the 
scattering matrix network model (Chalker Coddington model). 



B. Local Density Of States As Order Parameter 

Having established the role of r/L^ as the relevant scaling parameter for correlations of 
eigenstate amplitudes with universal exponent z(q) (related to r(q) by a scaling relation) let 
us now discuss the consequences of this for the interpretation of the local density of states 
being an order parameter of the LD transition. With the smearing-out of the 5-functions 
the LDOS is given by (see Eq. (M)) 



p(E,r) = r(E)- 1 \^(E,r)\ 2 (286) 

where \ip(E, r)| 2 stands for the microcanonical average of squared amplitudes at a given 
energy E. Since T(E) behaves as L~ d , the scaling behavior of the LDOS is determined by 
that of the wave function. Consequently, we have 

([p(E,r))«) L cxL^ d -^ (287) 

and for the typical value 

p t = exp[(ln(p(r))) J oc L d ~ a ° (288) 

which does reflect the LD transition. Scaling L with the correlation length £ c oc r~ v the 
typical LDOS vanishes on approaching the critical point with exponent /3 t = v{a.Q — d) (w 0.7 
in quantum Hall systems). 

The behavior of the typical LDOS is, according to the argument presented in the begin- 



ning of this section, similar to that of the participation ratio introduced by Wegner ||133 |. 
The participation ratio was based on the second moment of local probability. Higher mo- 
ments will work as well and our notion of typical LDOS is just a bit more convenient since 
it focuses on a typical value rather than on some specific moments. 

Our interpretation is also consistent with the findings for the LDOS in the localized and 



delocalized phase as discussed in Sect. VI. Therefore, we summarize the findings 
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1. In the localized phase the typical local density of states vanishes. 



2. In the delocalized phase the typical local density of states is finite and positive. 

3. Approaching the LD transition from the delocalized phase the typical LDOS vanishes 
with a positive critical exponent, 

A = K«o - d) . (289) 

A schematic picture for the typical LDOS that stresses its role as order parameter is shown 
in Fig. 31. 

The unconventional feature as compared to ordinary critical phenomena lies in the facts 
that 

1. The order parameter field has a broad distribution resulting in a non-linear dependence 
of exponents on the degree of moments considered (multifractality). 

2. The average value shows a vanishing scaling exponent while the typical value gives rise 
to a positive scaling exponent 

The scaling relations that we derived for the wave functions amplitudes transform to 
scaling relations of the LDOS since each box amplitude has to be multiplied by a constant 
factor of L d , 

((p(E, vi)) q (p(E + lu, p 2 ))*> £ oc (r/L^-^L-^ (290) 

with 

z(q) = d + 2r(q) - r{2q) , %) = 2(1 - q)d + 2r(q) . (291) 



These scaling relations form, in the sense of Sect. |V11 A| , the appropriate scaling relations 



connecting the spatial correlations of the local order parameter field to its scaling dimensions 
(cf. Eqs. (m,(|9D,(|9lD). 



We mention that i] := z(l) = d — D(2) « 0.4 ^ fj := z(l) = (for the correlator of 
the density of states in quantum Hall systems) with L^/r forming the scaling parameter 
is equivalent (cf. p8]J)7]]) to the phenomenon of anomalous diffusion found by Chalker and 
Daniell ||27|| . As pointed out in ||27|| , the anomalous character of diffusion lies in the non- 
Gaussian dispersion of a wave packet in time t despite the fact that the average diameter 
grows like This non-Gaussian time dispersion is caused by the multifractal character of 
eigenstates. 

A further support for the order parameter interpretation comes from conformal mapping 
ideas for 2D systems. In ordinary critical phenomena theory, the following assumption for 
critical correlation functions x{ r ) r ~ fl seems plausible: Scale invariance, reflected by power 
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laws, should hold also for local scale transformations which preserve angles but may change 
scales locally (called conformal mappings) 1221 . 



Thus, the main idea behind conformal mapping arguments in scaling theory is an ex- 
tension of a homogeneity law for correlation functions with respect to rescaling. Such a 
law exists also in the multifractal case, since any rescaling of all length scales in the LDOS 
correlator (for a fixed energy) by the same scaling factor s leads to 

(pi(sr)pi(sr')) sL = s~ s ^ (p%v)pi{v')) L . (292) 



where z(q) is given by Eq. ( |291 ). 



Extending this law to conformal mappings of a geometry Q to geometry Q leads, for 
large but finite 2D systems, to 

=M2l)l ' 111(22,1 ' • (293) 

where w(z) is any holomorphic function of complex coordinate z and w'(z) denotes the 
derivative. 

By choosing w(z) = (Lt/2tt) In z which maps the plane onto a strip one can show (cf. p2| ) 



that the correlator in the strip is characterized by a g-dependent quasi-lD localization length 
£(q; L t ) which is related to z(q) by 

£(q;L t )/L t = 2/(nz(q)). (294) 

In FSS calculations it is rather the average of the logarithm of the correlation func- 
tion which is calculated by knowing that this quantity defines a typical localization length. 



Therefore, the result of Eq. ( 294 ) leads in particular (consider q close to 0) to the relation 



A*= . 2 , d = 2. (295) 

Here A* appears since it is just the typical localization length at criticality divided by the 
width L t . Equation ( |294| ) has been confirmed in numerical calculations by Dohmen et al. 
], and Eq. (|295|) is, so far, in accordance with all known numerical results. 



C. Conductance As Scaling Variable 

The fact that the critical value «o — d is related to the critical value of the scaling variable 
(at least in 2D) is an important hint for the typical conductance g t being a one-parameter 
scaling variable. The reasoning behind this statement goes as follows: as pointed out in 
Sect. |VIII C| the scaling variable A(L t ) of finite size scaling in quasi-lD can be related to 
the typical conductance of a cube g t (L = L t ). From the 2 + e field theoretic calculations 
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at the LD transition ||133|| as well as from the multifractal LDOS found in 2D [|5(J we can 
conclude: if the correlation length is larger than system size and the typical conductance is 
still large, the eigenstates behave critical and are multifractal. The multifractal exponent 
«o — d is determined by g t and 



a - d oc l/g t . 



(296) 



At a generic critical point with A* ~ g t ~ 0(1) the conformal mapping relation tells, that 
also here the value of «o — d is determined by g t . This means that the bulk of the LDOS 
distribution function is, in the critical regime, already determined by the value of typical 
conductance. 

Before we collect results for the conductance distribution at criticality we discuss a similar 
object related to the level statistics. Not only the LDOS and conductance distributions 
should be universal functions at criticality (see Eq. ( |181| )), but also e.g. the level spacing 
distribution P(s; L) and the two-level correlation function R(s; L). To study these quantities 
in view of the one-parameter scaling picture was put forward in a work by Shklovskii et al. 



119|| . These authors pointed out that the study of the level spacing distribution P(s; L) is, 



starting from a Hamiltonian modeling, numerically much easier to perform than calculating 
conductances. Only the spectrum has to be calculated. According to the one-parameter 
scaling picture P(s; L) should obey the following scaling behavior 



lim P(s; L) 

L— >oo 



lim F(s; ocl) 

L—>oo 



Pi 



L,D,C 



(297) 



Here F(s; oll) is a function that does no longer depend on microscopic details, but only on 
one scale dependent parameter a L . Depending on the initial value of az the system flows 
to either the localized fixed point (L), the delocalized fixed point (D) or it stays at the 
critical fixed point (C). In terms of the typical conductance, ql = gt(L), the fixed points 
are characterized by 



lim g t (L) 

L—fOO 



(L) 

9t (C) 
oo (D) 



(29£ 



The corresponding asymptotic level-spacing distributions -Pl,d,c( s ) are given by the Wigner 
surmise, Eq. 



(299) 



P D (s) = A f) s p exp [-Bps 2 
the Poisson distribution, Eq. ( |120| ), 

P L (s) = exp (s) , 



(300) 
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and by a, yet unknown, critical distribution Pc{s). 

Numerical studies by several authors (e.g. | |119| |, [ f4"9"| , [pSfl , [ |135| |, |3J]) have demonstrated 



that a size independent critical distribution Pc{s) distribution exists, once the parameter 
r triggering the LD transition is put to its critical value r = 0. It is also known that, 
qualitatively, the distribution resembles the Wigner surmise for values s < so ~ 2 and the 
Poisson law for s > Sq. At So the distributions intersect. Furthermore, the way the critical 
distribution is approached allows for a calculation of the critical exponent of the correlation 
length v. This can be done in several ways. The authors of ||119|| took as a scaling variable 



the integral over P(s; V) up to s (one could use also some other point) 



so 

a(L;r) := J ds P{s; L;r) . (301) 



and analyzed it according to the general ansatz, (Eq. fll60|) , 

a(L; r) = a* + AtL 1/v + O (r 2 ) . (302) 



The values for v found by this method do agree with those found in FSS analysis within the 
errors. 

As to the precise form of Pc( s ) there exists no conclusive theory, except that the very 
small s behavior is again dictated by the level repulsion oc s 13 . For example, the tail of Pc( s ) 
is not exactly known. It seems now to be established that there is a leading exponential tail 
with sub-leading stretched exponential tail contributions. The question about the form of 
the tail is related || to the question about the relation between the number variance S 2 and 
the average level number (n), (see Sect. [VTD . Based on one-parameter scaling arguments 
Aronov et al. |fnj had predicted (for the LD transition point without magnetic fields) the 



following scaling law, for large (n), 

S 2 oc {{n)f~^ . (303) 

However, in this work the anomalous diffusion at criticality (or in other words the multifrac- 
tality) was not taken into account. Correcting this omission leads to an additional linear 
term, characterized by a finite level compressibility |33[ 



rfS 2 d-D(2) , , 

TTT = 7^ • 304 

d(n) 2d V ; 



In view of the finite level compressibility Eq. ( |303| ) describes sub-leading contributions to 
S 2 . 

It is very interesting that the analysis of S 2 as a function of (n) seems to allow for 
computing the correlation length exponent v and the fractal dimension D(2) related to the 
order parameter exponent (3 t = (ao — d) v (within the parabolic approximation to f{a) 
D{2) «d-2(ar o -<0)- 
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As to the critical conductance distribution the general considerations formulated in 
Eqs. (|297|,|298p f° r the level spacing distribution do apply as well. The precise form of 



the critical conductance distribution V*(g) is not known. Based on the 2 + e calculations by 
Al'tshuler et al. ||, Shapiro and coworkers [|36| predicted that power law tails should occur. 
So far, this could not neither be proved nor disproved by numerical means. However, in a 
work by Markos and Kramer |JT| it was shown that the critical conductance distribution is 
universal for a 3D system without magnetic field, i.e. it does not depend on system size L 
and does not depend on which parameter r is taken to its critical point value. The distri- 



bution showed fluctuations J ((Sg) 2 ) that are of the same order (0(1)) as the mean value of 
the distribution. The precise form of the tails could not be determined. 

The scaling of the conductance distribution close to criticality follows from expanding 
the scaling variable around criticality which brings the exponent v into play, 

V(g; L; r) = F(g; a L (r)) = V*(g) + F'{g- a*)ArL^ + ..., (305) 

where F' denotes the derivative of the function F with respect to the scaling variable ai- To 
analyze the behavior of V(g; L; r) at one particular point is numerically not advantageous. 
Therefore some integrated quantities work better. However, as explained earlier, the chosen 
quantity should not be sensible to the tails of the distribution. Therefore, one should refrain 
from taking moments. Unfortunately, not much has been done in this direction. 

Fastenrath et al. have investigated a quantity which is similar to the dissipative 
conductance, the so-called Thouless number, for quantum Hall systems and studied the dis- 
tribution function. Also in this case it turned out that, at the critical point, the distribution 
function was scale independent with fluctuations of the same order as the typical value. 
They explicitly took the typical Thouless number g t = exp (In g) and analyzed it according 
to 

g t {L;r)=gl + ArL 1 l v + ... . (306) 

They found a value for u 2.3) which is in good agreement with the value obtained by FSS. 
Thus, their observation gives further support that the typical conductance is an appropriate 
scaling variable. 

Cho and Fisher have performed a calculation of a true two-probe conductance within 
the Chalker-Coddington network model for the quantum Hall effect. There findings are 



consistent with the two-probe experiment by Cobden and Kogan |35j mentioned in the 
introduction: the distribution is almost uniform between and 1. This means in particular 
that the fluctuations are of the same order as the average value. For this uniform distribution 
on a finite interval the average value is already a convenient choice for the scaling variable 
and Cho and Fisher were able to show that the average value shows scaling of the form 
written in Eq. ( |306| ) with the correct critical exponent v « 2.3. Similar results were also 



obtained by Wang et al. [ |128| . 
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Thus, although the precise form of critical conductance distributions is not known, there 
are a number of results which confirm that it is universal and that the notion of typical 
conductance allows for studying the critical exponent of the correlation length. 

X. SUMMARY 

In this review we have discussed a few fundamental experiments that show the impor- 
tance of quantum mechanical interference effects for transport measurements in mesoscopic 
systems. Static disorder leads to complicated interference patterns of electronic wave func- 
tions which can cause large fluctuations in physical quantities. In addition, they can lead to 
localization of states to a finite volume within the mesoscopic system. Mesoscopic systems 
can be experimentally realized at very low temperatures (< IK) and with small devices 
(< 1/zm) such that the phase coherence length is larger than system sizes L. 

We introduced the basic physical quantities which are the local and global density of 
states (LDOS and DOS), the global conductance which is (in atomic units e 2 /h) given by 
a dimensionless quantity g, and the conductivity tensor a^ u (r,r';uj) (or equivalently the 
diffusion function D^ u {r, r';u;)), related to local charge transport. The relevant physical 
energy scales are the quantum kinematic scales of level spacing A and Fermi energy E F , and 
the transport energy scale set by the Thouless energy E Th . The conductance is given by the 
Thouless formula, g = E Th / A. A quantum kinematic length scales is the Fermi wavelength 
\f, and transport length scales are the microscopic mean free path l e and the localization 
length £ for localized states. 

We discussed models for disordered mesoscopic electron systems in terms of Hamilto- 
nians and scattering S'-matrices. From the Hamiltonian models energy eigenvalues e a and 
eigenvectors vfj a can be calculated. Their statistical properties determine the statistics of the 
most important quantity, the Green's function which determines all of the relevant quanti- 
ties and scales. Within the ^-matrix modeling the transmission strengths T a p are the most 
important quantities. They directly determine the conductance via the Biittiker formula, 
g = YlapTap. In addition, certain network ^-matrix models also allow for a determination 
of energy eigenvalues and eigenvectors. 

For the statistical and scaling properties of mesoscopic electron systems we collected the 
following insights. 

1. Ideal localized systems (with vanishing localization length £, i.e. vanishing conduc- 
tance g — 0) are characterized by a compressible spectrum of uncorrelated energies 
and a LDOS that typically vanishes, although the average DOS is finite. 

2. Ideal delocalized systems (with infinite Thouless energy E Th , i.e. infinite conductance 
g = oo) are characterized by an incompressible spectrum of correlated energies and 
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Gaussian distributed wavefunctions which leads to an exponential tail for the LDOS 
probability distribution. 

3. Changing control parameters, r, of the mesoscopic system like particle density, disor- 
der strength, pressure or applied fields can drive the system through a disorder induced 
localization-delocalization (LD) transition (at r = 0) which resembles a critical phe- 
nomenon. Due to the randomness of disorder the transition has to be described in 
terms of distribution functions the flow of which with respect to system size (scaling) 
forms the object of interest in a critical phenomenon description. 

4. A one-parameter scaling theory for distribution functions of global quantities like con- 
ductance and energy-level spacing seems to work well. It states that the distributions 
are essentially determined by one scale dependent parameter ai, called scaling vari- 
able, that determines the flow, V(X; L) xs F(X; a L ). 

5. A possible scaling variable for the LD transition is given by the geometric mean of 
the conductance distribution which we called typical conductance g t . With increasing 
system size g t can approach the idealized localized phase (g t = 0), or it can approach 
the idealized delocalized phase (gt = oo). Corrections to distribution functions as 
compared to the ideal situations are controlled by finite values of g t . Right at the LD 
transition g t will become a size independent universal quantity gl that is of order 1 in 
generic LD transitions. The corresponding distribution functions of global quantities 
like conductance and level-spacing are universal scale-independent functions. 

6. The fictitious system size £ c for which <7t(£c) ~ 9* defines the correlation length of the 
LD transition. In the localized phase it can be identified with the localization length 
£. Close to the transition point, where it diverges, it scales with a critical exponent, 

OC T-". 

7. While in the localized phase the typical (not the average) LDOS vanishes, the LDOS 
shows power law scaling in non-localized phases. The average value is always scale- 
independent. In the ideal delocalized phase one exponent, simply related to dimension- 
ality d, controls the scaling of the whole distribution function. At the LD transition a 
spectrum of scaling exponents is necessary to describe the scaling of the distribution 
which is close to a log-normal distribution. The spectrum has an interpretation in 
terms of a spectrum of fractal dimensions (multifractal spectrum). Still, there is one 
particular exponent, related to the typical value of the LDOS, which determines the 
log-normal approximation for the distribution. It is called — d > and is deter- 
mined by the value of the scaling variable at criticality. Thus, multifractality is not in 
conflict with one-parameter scaling theory. 
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8. On the contrary, the multifractal scaling of the LDOS opens the possibility to consider 
the typical value of the LDOS as an appropriate order parameter for the LD transition. 
The corresponding critical exponent is given by [3 t = v{oiq — d). Further support to 
this idea comes from the observation that a number of scaling relations for the LDOS 
can be formulated that are in close analogy to similar relations in ordinary critical 
phenomena. 

In short, the occurrence of broad distributions and mult ifr act ality in disordered meso- 
scopic electron systems is by no means in contradiction to the one-parameter scaling theory, 
but points out that only typical rather than average quantities can serve as order parameter 
and scaling variable in this theory. 
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APPENDIX A: BRIEF ACCOUNT OF LINEAR RESPONSE THEORY 

We start with a general outline of how to calculate linear response quantities when 
the system is specified by a Hamiltonian. We then specialize to additive systems where 
the resolvent G(z) := (z — H)' 1 of the corresponding one-particle Hamiltonian (taken at 
complex energies z) yields all the basic physical quantities. Studying transport in multi- 
probe geometries will lead to the concept of the scattering matrix. 

Consider a many-particle system characterized by the Hamilton operator H coupled to 
a linear increment of a force field, F, 

Tip — H — XF (Al) 

where X is a hermitean operator of the system. The system is assumed to be in an equilib- 
rium state at temperature T described by the (grand) canonical density operator p F . The 
thermodynamic susceptibility x Y x * s defined as 

Linearizing pp with respect to the increment F the thermodynamic susceptibility reads 
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Xyx =P(AX;AY). (A3) 
where the bilinear expression (Kubo scalar product) 

(X; Y) := - J dX (XY{ih\)) = (Y; X>* (A4) 



/5 

u 

is called the canonical correlator of X and Y with /? = (fc^T) -1 , AX := X — (X), and X(t) = 
e (t/h)Ht-^ e -(i/h)Ht - g ^j me evolution of X in the Heisenberg picture. The isothermal 
susceptibility describes an idealized linear response situation where the force field does not 
disturb equilibrium, but is just a parameter field of the system. 

We turn now to the dynamic linear response problem @,[5(J where the following (Kubo) 
process is considered: the increment in an applied force field F(t) depends on time t and is 
switched on adiabatically at t = — oo. The total Hamilton operator is 

n F (t) = H - XF(t)e"* , (A5) 

where the slowness parameter 77 has to be much smaller than the inverse of a typical positive 
time scale t m at which the response to the field is measured The system is assumed to be 
in equilibrium for t = —00 described by the density operator p(F = 0) and isolated for 
t > —00. Restricting our consideration to one Fourier component of the force field, 

F(t) = Fe~ luJt , (A6) 

the resulting linear increment in the Fourier component of the mean value of an observable 
Y is 

(Y) F = x YX (z)F (A7) 
where the dynamic susceptibility x Y x( z ) ^ s gi ven by the Kubo formula 

Xyx (z) = P (X II Y) [z] = fi (AX II AF) [z] . (A8) 
Here X denotes the time derivative of X and the dynamic correlator 



(X II Y) [z] := J dte izt (X;Y(t)} (A9) 


appears as the Laplace transform (z = u + irf) of the time dependent canonical correlator. 

Integrating by parts in Eq. (|A8|) yields the fundamental equation of motion for dynamic 
correlators 

(X || Y) [z] = (X; Y) + iz (X || Y) [z\ = - (X || Y) [z] . (AlO) 



86 



Another important equation of motion, referred to as Kubo relation, holds for the canonical 
correlator 

(X;Y) = ^([X,Y\) (All) 

where [X, Y] denotes the commutator of X and Y. It is worth mentioning that the second 
term in Eq. (|A10| ) can be finite in the case of zero frequencies, u = 0, since it defines the 



long time average of (X; Y(t)). 



(X; Y) := lim rj (X || Y) [irj] . (A12) 



Equations (|A"3,X8, AlC , All ) form the backbone of any linear response theory (cf. 



With these we can now write down the expressions for the physical quantities considered in 
Sect. [n|. 

The local charge density, q(r), can be related to local electric potentials, U(r), by the 
charge response function, Il(r, r u)Q 

q(r) = J d d r'U(r,r';u)U(r'). (A13) 

Charge response and conductivity (Eq. (|10D ) are related due to the continuity equation 

divj(r) -iuq(r) = 0, (A14) 
d d 

S^^^ 1- '^ =iwn(r,r';w). (A15) 

fj,v u ' A* u ' v 

The particle density n(r) couples to the local chemical potential /x(r) in the same way 
as the local charge q(r) = —en(r) couples to the local electric potential U(r) and thus the 
thermodynamic charge response reads 

n(r,r') = -/9(Ag(r');Ag(r)), J d d r' U{r, r ') = e 2 p{r; p) (A16) 

while the dynamic charge response reads 

n(r, r '; uj) = -13 (Ag(r ') || Ag(r)) [z] = U(r, r ') - ij3z (Ag(r ') || Ag(r)) [z] . (A17) 

Note, that the equilibrium charge response Il(r,r') is not equal to the uo — ► limit 
of the dynamical charge response II(r, r';u). This can be also seen when considering the 
response to homogeneous potentials. While a homogeneous electric potential does not create 
a charge response, and hence / d d r'U(r, r';u) = 0, shifting the overall chemical potential by 
a constant usually leads to a change in the particle density, i.e. / d d r' Il(r, r') = e 2 p(r; fi). 



9 For simplicity we drop the dependence on the time dependent vector potential. 
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According to the general scheme described above, The conductivity tensor appears as 
the dynamic current correlator 



a tlu (r,r';u;) = P(j u (r>) \\j„(r))[z] (A18) 

and the relation in Eq. ( |A14 ) can be viewed as a special case of Eq. (|A10[) . The corresponding 
global conductivity tensor reads in the thermodynamic limit L — > oo 

e 2 

a^(to) = —[3{v u \\v^[z\ (A19) 

where v denotes the velocity operator of electrons. 

More delicate is the representation of the global conductance coefficients Gk m (uj), since 
one has to know the local potential distribution within the conductor as a functional of 
applied electro- chemical potentials U m and we do not go into the detail of this interesting 
problem. However, due to the continuity equation the situation simplifies for the d.c. situa- 
tion: The current in any lead is independent of the actual potential distribution U(r) within 
the conductor (cf. ]7D| ) and Gkm is determined by the dynamical correlator of asymptotic 
current operators which are the time derivatives of the total charge operator, Ik = Qk, 

G km = -(3 (Q m || Q k ) (A20) 

where omission of the argument [z] refers to the d.c. situation. 

Still, the calculation of the correlators requires the knowledge of the many particle density 
operator Pf=o and the machinery of many-body physics is necessary to get quantitative 
results. 

Assuming additive electron Hamilton operators like the ones introduced in Sect. [IV] the 
correlators can be written in closed form with the help of the one-particle resolvent (Green's 
function in real space representation) 

G ± (E) := (E±ie-H)- 1 , e -> +0 (A21) 

of the corresponding one-particle Hamiltonian H. It defines a spectral resolution via 

5{E-H) = {2m)- 1 {p-{E) - G + (E)) . (A22) 

With the help of the Fermi distribution function 

f(E) = [l + e^-")]" 1 , (A23) 

(where [i denotes the chemical potential) average values can be calculated by 

(X) = ^ J dE f(E)Tr {G-(E)X - G + (E)x} . (A24) 
The canonical correlator reads 
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(AX; AF) = 7^jJ dE f(E)Tr \XG~YG~ - XG + YG + ) (E) , (A25) 
and identifying 77 with 2eh~ 1 yields 
(X || Y) [z) = ^J dEf(E)Tr{(G- - G + ) e (XG^G~ E Y - Y G E+hu) G+x)} • (A26) 

for the dynamical correlator, provided the equilibrium mean values of X or Y vanish. 

Since the Fermi function is a step function at T = the general structure / dE f(E)A(E) 
occurring in Eqs. ( |A24HA26|) allows to conclude e.g. for the canonical correlator 

df 



(A27) 

T=0,fi=E F =E 



(3(X;Y) =-J dE^(3(X;Y) 

A similar equation holds for the dynamic correlator. Therefore, for additive (non-interacting) 
Fermion systems we can always restrict the discussion to the zero temperature limit. 

In the zero frequency limit further simplifications are possible. For example, any auto- 
correlator is essentially determined by states at the Fermi level E F , 

(3 (X || X) (T = 0) = tiTi Tr {X5(E F - H)X5(E F - H)} (A28) 

which shows that the longitudinal conductivity is a Fermi level quantity. 

To simplify other correlators it is convenient to use the equation of motion on the level 
of the one-particle resolvent operator 

G(z)XG(z') = % - {G{z)X - XG(z') + {z - z') G(z)XG{z')} . (A29) 
n 

An important example for an (essentially) Fermi level quantity is the zero frequency dynamic 
correlator of two operators one of which is a time derivative 



(3 (X || Y) (T = 0) = (3(AX; AY)(T = 0) + lim - Tr \XG-(E F )YG + (E F )\ . (A30) 

v ' e— >+0 7f I. > 



The d.c. conductance coefficients at T = read for n 7^ m 

G km {u = 0,T = 0) = ^ Tr {G-(E F )Q m G + (E F )Q n } (A31) 

and it is obvious from Eq. (|A30| ) that the conductance coefficients are Fermi level quantities, 
too. 

By applying Eq. ( |A30| ) to the longitudinal conductivity one finds 

a xx {u = 0, T = 0) = -j^e 2 Tr {tTG+x 2 - CTa;G + x} (E F ) . (A32) 

From Eq. (|A31|) one can conclude that the conductance coefficients are determined by 
asymptotic Green's functions, and it is tempting to formulate a scattering theoretical de- 
scription. This is indeed possible by decomposing the total Hamiltonian into two parts, one 
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of which describes the leads, Hi eads , and a second part which describes the conductor plus 
its coupling to the leads, H scatter , 

H = #lcads + -^scatter • (A33) 

The total resolvent G(E) can then be decomposed according to Dyson's formula on intro- 
ducing the T-matrix of scattering theory 

G(E) = G°(E) + G°(E)T(E)G°(E) , T(E) = H scaittcr G(E)(G°(E)y 1 . (A34) 

The eigenstates of the leads are scattering states, denoted as \k,a), where k labels the lead 
and a the quantum number (channel) of a scattering state at energy E. The matrix elements 

t a k i(E):=(m,P\T(E)\k,a) (A35) 

define transmission amplitudes to scatter from channel a in lead k to channel (3 in lead m 
and form the matrix of transmission amplitudes t^ n . Applying Eq. ( |A30|) and Eq. ( A31 ) and 
using the fact that the asymptotic charge operators are nothing but — e times the projector 
onto the Hilbert space of the leads, the important Buttiker formula (m 7^ k) pOj 

e 2 

G km (cu = 0, T = 0) = - Tr { WL) ( e f) ( A 36) 

can be obtained. 
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Figure Captions 



Figure 1: Schematic view of a mesoscopic conductor. Electrons experience elastic scattering 
events on their way from the current source to the current sink. 

Figure 2: A mesoscopic ring-shaped conductor. A flux is put through the ring. 

Figure 3: Resistance oscillations for a ring-shaped conductor with flux 0. The period is 
(f)o = h/e as shown in the inset where the Fourier spectrum is displayed, (after ||124j| ). 



Figure 4: Universal conductance fluctuations as a function of applied magnetic field. Differ- 
ent curves correspond to gold-nanowires of different length. The average value of conduc- 
tance was shifted for convinience (after 

Figure 5: Hopping transport in amorphous Silicon. The logarithm of the conductivity is 
plotted versus T _1//4 showing Mott's law. The different curves are for different deposition 
temperatures (after [pp. 

Figure 6: The resistance of thin disordered film of coupled fine Cu-particles as function of 
logarithm of the temperature (after 



Figure 7: Diffusion path in a disordered system. The electron propagates in both directions 
(time reversal symmetry). Due to constructive interference the quantum return probability 
is twice as great as the classical return probability. 

Figure 8: The resistance curves of a thin Mg-film (upper set of curves) as function of 
magnetic field (magneto- resistance). The weak localization occurs for zero field. After a 
superposition with 1/100 atomic layer of Au the magneto-resistance changes drastically. 
The Au introduces a rather pronounced amount of spin-orbit scattering which rotates the 
spin in the complementary scatterd waves. This changes the interference from constructive 
to destructive (after [|T8|1 ). 

Figure 9: Conductivity as a function of electron concentration indicating strong localization 
(after fl22j ). 

Figure 10: The quantum Hall effect: The Hall resistance p xy and the magnetoresistance p xx 
as function of the magnetic field measured in medium mobility GaAs heterostructures at 
60mK (after |2j] ). 

Figure 11: Schematic view of a multi-probe conductor. Currents Ik entering the system 
from probe k are counted positive. 
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Figure 12: Schematic view of a box shaped conductor in d dimensions characterized by 
length L and cross section Lf~ x . 

Figure 13: Schematic view of energy branches for semi- infinite leads. Due to transversal 
quantization a discrete number of modes exist (characterized by a wave number k) for a 
given value of the Fermi energy Ep. 

Figure 14: Scattering in a quasi-lD system. 

Figure 15: Adding two ID conductors in series. 

Figure 16: The scattering element of a 2D network. For unit amplitude on the ingoing link 
from the left and all other ingoing links being empty the amplitudes on the outgoing links 
are r\L, tr and r, respectively. 

Figure 17: The network resulting from the elements of Fig. 16. 

Figure 18: The basic scattering elements of the Chalker-Coddington model. 

Figure 19: The Chalker-Coddington network model for 2D disordered electrons in strong 
magnetic fields. For strong reflection to the left (right) the electrons perform chiral loops. 

Figure 20: Qualitative behavior of the /3-function for disordered electrons. For d = 2 the 
full line corresponds to absence of magnetic fields and spin-orbit scattering while the dotted 
line corresponds to the presence of spin-orbit scattering. In d = 3 a LD transition occurs 
quite generally at some critical value g* which forms an unstable fixed point of the flow 
determined by the /3-function. 

Figure 21: The /3-function as following from Eq. ( |169| ). The curves correspond to (from 
below) d = 1,2,2.2,3. 

Figure 22: Qualitative behavior of multifractal exponents. The top figure shows the expo- 
nents r(q) corresponding to moments of a multifractal distribution, the middle shows the 
resulting generalized dimensions D q = r(q) /(q — 1) and the bottom shows the /(a)-spectrum 
being the Legendre transform of r(q). f(a) describes the scaling of the whole distribution. 

Figure 23: Qualitative behavior of the conductance distribution in a metal. Within a range 
of width d ~ J (g) around the average value (g) the distribution is Gaussian with universal 
variance of order 1. Outside this region the distribution develops log- normal tails. 

Figure 24: The scaling function for a 3D quantum Hall system showing a LD transition. 
A(E,M) denotes the scaling variable as a function of energy E and system size M. £(E) 
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denotes the correlation length obtained from a fitting procedure. Different data symbols 
correspond to different energies and system sizes (after PDJ). 



Figure 25: Qualitative picture visualizing the explanation of the quantum Hall effect being 
due to a sequence of localization-delocalization transitions occurring at the Landau energies. 
p, <r xx , a yx denote the density of states, the dissipative conductivity and the Hall conductivity, 
respectively. In finite systems the range of extended states (grey) on the energy scale has 
finite width. This width is believed to shrink to zero in the thermodynamic limit. 

Figure 26: Squared amplitudes of a critical wave function for a system of linear size 200 
magnetic lengths. The corresponding f(a) spectrum (•) together with the parabolic ap- 
proximation (• ■ ■) and a semi-elliptic approximation (-) are also shown (after ||110| ). 



Figure 27: Histogram of the logarithm of squared amplitudes shown in Fig. 26. The con- 
tinuous curve is the distribution function following from the corresponding f(a) spectrum 
with a = 2.28 ± 0.02 (after PD| ). 



Figure 28: Squared amplitudes of a critical wave function in the Chalker-Coddington network 
model. (after |jT"5|| ). 

Figure 29: Histogram of the logarithm of squared amplitudes shown in Fig. 28. The con- 
tinuous curve is the distribution function following from the parabolic approximation to 
the f{ot) spectrum with a = 2.27 ± 0.02. The inset shows the histogram for the squared 
amplitudes itself (after 



Figure 30: The function of critical exponents z(q) following from the the scaling relation 
Eq. (|27|) (after jTTOlQ . 

Figure 31: Qualtative behavior of the typical local density of states reflecting a localization- 
delocalization transition at critical Fermi energy E*. 



99 




1 




2 




3 




4 



10-2 



o/(ficm )- 1 




I - - ; - ■■ | I 

0.22 0.24 0.26 0.28 

(K/T)"* 



5 




6 




7 




8 




9 



3000 




BITesla] 



10 




11 



E A 




13 



I 




I 1 




> 




< 





s 


o 1 


< 




> 



14 




16 



17 



T 

I 




19 




20 




22 




23 



> 



0.8 



t> 



*- 



0.4 - 



0.0 



-0.4 



-0.8 



O 



7* V 



s> o o 



-1.2 



-1.2 -0.9 -0.6 -0.3 0.0 0.3 0.6 0.9 1.2 



In (£(E)/M) 



24 



n+1 



A a xx °yx [e 2 /h] 



n-1 



A ft A 



A A 




25 




26 




27 




28 



N/N total 




-22 -20 -18 -16 -14 -12 -10 -8 -6 

InP 



29 




30 




31 



